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A  theory  is  presented  for  the  exact  synthesis  of  five-,  six-  and 
higher-link,  planar  mechanisms  for  five  multiply  separated  position 
dynamic  specifications  (the  class  of  prescribed  conditions  for  synthe- 
sis consisting  of  inertia  forces  and  torques,  momenta  or  combinations 
of  positions  and  time  derivatives  up  to  the  fourth  order,  i.e.,  velocity, 
acceleration,  jerk  and  jounce) .   Also,  certain  previously  unsolved 
geometric  synthesis  problems  are  solved  utilizing  this  theory. 

A  generalized  computer  program,  MECSYN,  is  presented  for  the 
synthesis  of  multi-link,  planar  mechanisms  to  allow  practical  use  of  the 
theory  above  in  addition  to  the  usual  synthesis  for  geometric  specifica- 
tions.  MECSYN  is  a  closed  form,  noniterative  program  which  uses  a 
modular  approach  employing  four  tools  of  synthesis  (constraint  link 
synthesis,  inversion,  angular  cognates  and  path  cognates)   for  all  pos- 
sible cases  of  four  and  five  multiply  separated  positions.   Other  fea- 
tures of  MECSYN  are  complete  mechanism  analysis  for  the  specified  pre- 
cision conditions,  link  length  ratios,  Grashof 's  criteria,  straight  line 
motion  synthesis,  and  function  generation  for  four-,  five-  and  six-link 
mechanisms . 

ix 


Examples  are  presented  for  four-,  geared  five-  and  six-link 
mechanisms,  including  two  practical  examples.   A  Stephenson  3  and  a 
Stephenson  2  six-link  mechanism  is  designed  to  transfer  an  object  from 
a  continuously  moving  conveyor  to  a  continuously  rotating  turret  so 
that  impact  forces  are  eliminated.  A  Stephenson  3  mechanism  is  designed 
to  allow  automatic  torch  or  bandsaw  cutting  of  large,  low  pitch 
(Pd  =  3/8)  gear  teeth. 


CHAPTER  I 


INTRODUCTION 


A.  A  Short  History  of 

"Kinematic"  Synthesis 

In  recent  years,  many  researchers  have  investigated  the  problem 
of  approximating  a  function  or  the  coplanar  motion  of  a  body  with 
a  mechanism.   The  impetus  for  this  research  was  the  advent  of  the  digital 
computer.   Although  some  of  these  solutions  had  been  obtained  by  graph- 
ical means  in  previous  years,  many  required  enormous  efforts  and  pains- 
taking attention  to  accuracy  to  produce  mechanisms  for  the  maximum 
number  of  solvable  specified  conditions.   Also,  generation  of  a  large 
number  of  general  mechanisms  to  allow  selection  on  the  basis  of  pre- 
ferred mechanism  properties  was  usually  impractical. 

Most  of  this  research  was  concentrated  on  synthesizing  a  four- 
bar  mechanism  such  that  a  point  in  the  coupler  plane  passed  through 
three  to  five  positions  while  the  coupler  angle  had  prescribed  rotations 
between  these  positions.   This  is  referred  to  in  the  literature  as 
finitely  separated  positions  (FSP) ,  while  specification  of  geometric 
derivatives  for  the  path  point  is  referred  to  as  inf initesimally  sepa- 
rated positions  (ISP).   No  derivative  conditions  were  specified  for  the 
coupler  rotation  angle  since  it  was  the  independent  geometric  parameter 
for  the  specified  path  point  derivatives. 


This  type  of  synthesis  has  been  known  in  the  literature  as 
"kinematic"  synthesis  when,  in  fact,  it  has  not  included  specifications 
for  kinematics  in  the  fullest  definition  of  the  term  (such  as  velocity, 
acceleration,  jerk  and  jounce).   The  objective  was  to  synthesize  mechan- 
isms to  satisfy  purely  geometric  properties  (such  as  the  tangent,  curva- 
ture, inflection,  etc.  of  a  given  function).   The  result  of  this  type  of 
synthesis  is  a  mechanism  which  satisfies  certain  geometric  properties 
regardless  of  how  it  is  driven  (i.e.,  input  velocity,  acceleration, 
jerk,  etc.),  and  will  be  referred  to  hereafter  as  geometric  synthesis. 

Five  finitely  separated  position  synthesis  of  four-bars  was 
done  graphically  by  Burmester  [1]  in  1888.   This  problem  was  solved 
analytically  by  Primrose,  Freudenstein  and  Sandor  [2]  in  1964  and 
Bottema  [3]  also  in  1968.   Graphical  studies  on  five  infinitesimally 
separated  positions  were  given  by  Mueller  (1895,  [4]),  Allievi  (1895, 
[5]),  Krause  (1920,  [6]),  Wolford  (1960,  [7])  and  Veldkamp  (1963,  [8]). 
A  unified  theory  for  treating  all  seven  combinations  of  finitely  and 
infinitesimally  separated  positions  was  developed  by  Tesar  et  al.  (1968, 
[9]),  and  is  referred  to  as  five  multiply  separated  positions  (MSP). 

Frequently,  the  number  of  free  parameters  for  synthesis  avail- 
able in  four-bar  mechanisms  is  not  sufficient  to  adequately  solve  the 
problem  at  hand.   In  this  case,  the  next  two  classes  of  constrained 
planar  mechanisms  that  can  be  synthesized  are  geared  five-link  and  six- 
link  mechanisms.   Dynamic  specification  synthesis  as  presented  here 
would  not  be  possible  without  the  capability  to  synthesize  six-link 
mechanisms  due  to  the  large  number  of  parameters  required. 


Function  generation  for  the  Watt  2  (two  four-bars  in  series) 
has  been  investigated  several  times.   Solutions  have  been  presented 
for  five  finitely  separated  precision  points  in  closed  form  and  numer- 
ical techniques  have  been  presented  for  up  to  nine  precision  points  by 

authors  such  as  Kramer  and  Sandor  [10],  Levitskii  [11]  and  McLarnan 

* 
[12].   Hain  presents  several  synthesis  techniques  for  six-bar  linkages 

[13].   However,  all  the  solutions  are  obtained  by  graphical  techniques 
and  vary  from  two- to  five-position  synthesis. 

A  uniform  theory  for  the  solution  of  more  than  35  distinct  six- 
link  synthesis  problems  for  five  FSP  and  in  many  cases  five  MSP  was 
presented  by  the  author  [14]  and  the  author  and  Tesar  [15]  in  1971. 
Four  basic  tools  of  synthesis  were  utilized  (curvature  transform, 
inversion,  path  cognates  and  angular  cognates)  of  which  only  the  first 
two  had  been  used  extensively  in  the  past.   This  work  allows  the  solu- 
tion of  the  dynamic  specification  synthesis  problem  considered  in 
Chapter  II  by  providing  the  necessary  number  of  free  synthesis  parameters 
available  in  terms  of  six-link  mechanisms. 

Nine  precision  point  synthesis  of  geared  five-bar  mechanisms 
was  solved  numerically  by  Roth  [16]  in  1962  and  geared  five-bar  function 
generation  has  been  developed  by  Oleksa  and  Tesar  [17]  and  independently 


*It  should  be  noted  that  the  Watt  2  mechanism  has  no  coupler  motion  of 
higher  order  than  a  four-bar  (sixth  order  algebraic  curve),  while 
other  six-link  mechanisms  have  up  to  fourteenth  order  algebraic  curve 
coupler  paths.   As  a  result  the  Watt  2  is  the  simplest  to  synthesize 
and  has  restrictions  similar  to  a  four-bar.   One  such  restriction  is 
that  for  higher  than  five  precision  points,  the  existence  of  a  solu- 
tion is  not  guaranteed. 


by  Sandor  and  Erdman  [18]  in  1970.   The  author  solved  the  remaining 
problem  of  five-position  coplanar  synthesis  for  geared  five-bars  also 
in  [14]  and  with  Tesar  [19]  in  1971. 

All  of  the  preceding  references  solve  only  the  geometric 
synthesis  problem.   The  objective  here  is  to  synthesize  mechanisms  for 
dynamic  specifications  as  will  be  described  in  the  following  section. 

B.   Scope  of  This  Study 

Preliminary  work  in  the  synthesis  of  mechanisms  to  satisfy 
dynamic  specifications  was  done  by  Eschenbach  and  Tesar  [20,21]  and 
Eschenbach  [22]  for  a  four-bar  mechanism  by  a  different  method  than 
will  be  presented  here.   Also,  emphasis  was  placed  on  selection  of 
constraints  based  on  force  transmission  angle  of  a  four-bar.   Eschenbach 
[22]  concluded:   "Thus  it  is  impossible  to  coordinate  the  complete 
dynamic  state  of  a  simple  crank  and  body  beyond  the  second  order  accel- 
eration state,"  in  the  context  of  four-bar  mechanisms.   The  first  objec- 
tive of  this  study  will  be  to  give  closed  form  solutions  for  the  coordi- 
nation of  a  simple  crank  input  with  the  complete  dynamic  specification 
of  a  body  for  seven  possible  combinations  of  specifications  (described  in 
Chapter  II)  including  fourth  order  time  derivatives.   Because  of  the 
previously  noted  limitation  in  terms  of  available  parameters,  this  must 
be  done  by  synthesis  of  geared  five-link  and  six-link  mechanisms.   The 
objectives  may  be  stated  as  follows: 

1.   Presentation  of  a  theory  for  the  exact  synthesis  of  five-,  six- 
and  multi-link  mechanisms  for  dynamic  specifications  (Chap- 
ter II).   The  term  dynamic  specifications  is  used  to  mean  the 
class  of  prescribed  conditions  for  synthesis  consisting  of 


inertia  forces  and  torques,  momenta  or  combinations  of  posi- 
tion and  time  derivatives  up  to  fourth  order  (i.e.,  velocity, 
acceleration,  jerk  and  jounce).   This  does  not  include  all  pos- 
sible problem  definitions  where  a  general  applied  force  or  torque 
is  specified.   The  maximum  number  of  conditions  which  will 
allow  exact  solutions  are  specified  for  each  link  (i.e.,  five- 
position  problem).   Also,  certain  previously  unsolved  geometric 
synthesis  problems  are  now  solvable. 

2.  Development  of  a  generalized  computer  program  (Chapter  V)  for 
the  synthesis  of  four-,  five-,  six-  and  multi-link  mechanisms 
to  allow  practical  use  of  the  theory  listed  above  in  addition 
to  synthesis  for  purely  geometric  specifications.   A  modular 
approach  is  used  employing  the  four  tools  of  synthesis  employed 
by  the  author  in  six-link  synthesis  [14]  (constraint  link  syn- 
thesis, inversion,  angular  cognates  and  path  cognates).   Other 
synthesis  tools  useful  to  designers  are  included  such  as  analy- 
sis of  the  resulting  mechanisms  at  the  specified  precision  con- 
ditions, link  length  ratios  and  Grashof's  mobility  criteria. 

3.  Application  of  1  and  2  in  the  solution  of  examples  for  four-, 
five-  and  six-link  mechanisms  (Chapter  III)  including  two  prac- 
tical examples  (Chapter  IV). 

C.   Six-Link  Mechanisms 


A  constrained  closed  chain  is  one  in  which  the  motion  of  each 
link  is  uniquely  determined  with  respect  to  the  other  links.   This  study 
will  deal  only  with  constrained,  planar  mechanisms  having  revolute, 
prismatic  or  gear  pairs  (joints).   Most  mechanisms  can  serve  a  useful 


purpose  only  if  the  motion  is  predetermined,  given  a  single  input. 
On  occasion  mechanisms  are  desired  which  require  two  or  more  inputs  to 
produce  a  certain  motion  such  as  control  devices  or  differential  mechan- 
isms.  These  are  not  constrained  mechanisms  but  may  be  synthesized  by 
the  methods  presented  here. 

Due  to  the  work  of  Grubler  and  others,  it  is  possible  to  deter- 
mine the  number  of  degrees  of  freedom  of  a  planar  mechanism  from  the 
number  of  links  and  the  number  of  joints  alone, 

F  =   3(n-l)  -  2g  -  h, 

where  F  is  the  number  of  degrees  of  freedom,  n  is  the  number  of  links 
including  the  fixed  link,  g  is  the  number  of  lower  pairs  and  h  is  the 
number  of  higher  pairs.   Thus  for  a  closed  planar  turning  pair  chain 
to  be  constrained  (F = 1) , 

g=  (3/2)n  -  2, 

and  for  closed  six-link  turning  pair  mechanisms  the  number  of  turning 
pairs  must  be  seven. 

Grubler  also  proved  that  the  minimum  number  of  binary  links 
(having  two  kinematic  pairs)  in  any  constrained  chain  is  four  and  the 
greatest  number  of  pairs  by  which  any  link  of  the  chain  can  be  joined 
is  n/2  (n— total  number  of  links).   Thus  a  six-bar  chain  (constrained) 
must  have  at  least  four  binary  links,  cannot  have  any  quaternary  links 
(having  four  kinematic  pairs)  and  can  have  at  most  two  ternary  links. 

There  are  two  chains  having  four  binary  and  two  ternary  links. 
The  Stephenson  chain  Fig.  (1-1)  has  a  binary  link  separating  the  two 


Stephenson  3 
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Figure  1-1  Stephenson  Six-Link  Chain 


Figure  1-2  Watt  Six-Link  Chain 


ternary  links.   The  Watt  chain  Fig.  (1-2)  has  two  ternary  links  joined 
together.   There  is  one  chain  having  five  binary  links  and  one  ternary 
link  requiring  one  double  joint  (two  turning  pairs  having  a  common  axis 
of  revolution).   Six  binary  links  having  two  double  joints  make  up  the 
last  six-link,  constrained  chain.   The  last  two  chains  can  be  described 
as  special  cases  of  the  Stephenson  and  Watt  chains.   They  have  fewer 
free  parameters  for  synthesis  and  will  not  be  separately  considered  here. 

All  possible  distinct  inversions  of  the  Stephenson  and  Watt 
chains  are  illustrated  in  Figs.  (1-1,  1-2).   Some  inversions  of  a  chain 
are  equivalent  due  to  the  conceptual  symmetry  of  the  chain. 

There  are  three  distinct  inversions  of  the  Stephenson  chain. 
Fixing  binary  link  1  or  6,  Fig.  (1-lc)  produces  a  Stephenson  1  mechanism. 
Fixing  link  3  or  4  gives  a  Stephenson  2  mechanism,  Fig.  (1-lb) .   Fixing 
either  ternary  (2  or  5)  yields  a  Stephenson  3  mechanism,  Fig.  (1-la) . 

Fixing  any  binary  link  of  the  Watt  chain  produces  a  Watt  1 
mechanism,  Fig.  (l-2b) .   Fixing  either  ternary  link  gives  a  Watt  2 
mechanism,  Fig.  (l-2a)  which  is  equivalent  to  a  pair  of  four-bars  in 
series  (possessing  a  common  link  and  fixed  pivot  of  that  link) . 

D.   Four  Tools  of  Synthesis 

It  will  be  shown  that  the  objectives  stated  above  can  be 
achieved  by  the  use  of  four  tools  of  synthesis:   constraint  link  syn- 
thesis, inversion,  path  cognates  and  angular  cognates. 
1.   Constraint  Link  Synthesis 

Constraint  link  synthesis  has  also  been  termed  curvature  syn- 
thesis from  classical  curvature  theory.   This  is  probably  the  most 
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important  and  most  difficult  aspect  of  the  synthesis  process.   Given 
a  set  of  desired  specifications  for  a  link  or  links,  the  objective  is 
to  find  solutions  for  constraint  links  which  will  form  a  constrained 
mechanism  to  satisfy  the  specifications.   Several  different  methods  of 
constraint  link  synthesis  have  evolved  for  geometric  synthesis;  i.e., 
graphical,  analytical  and  iterative  numerical  techniques.   Although 
analytical  solutions  are  the  most  desirable  for  many  reasons  (accuracy, 
computing  time,  etc.)  they  are  the  most  difficult  to  formulate.   Since 
geometric  synthesis  is  not  sufficient  to  design  mechanisms  for  dynamic 
specifications,  a  new  curvature  transform  must  be  developed.   Chapter  II 
presents  a  closed  form  analytical  constraint  link  synthesis  formulation 
for  dynamic  specifications  also  called  the  dynamic  specifications  curva- 
ture transform  or  the  dynamic  curvature  transform.   Solutions  are  found 
for  the  maximum  number  of  specifications  in  terms  of  the  number  of 
parameters  for  the  mechanisms  treated. 

In  keeping  with  the  closed  form  analytical  nature  of  the  con- 
straint link  synthesis  presented  here,   all  other  synthesis  tools  are 
solved  in  closed  form  including  the  analysis  (Chapter  III).   This  allows 
the  computer  program,  MECSYN  (Chapter  V)  to  be  precisely  verified  and 
highly  efficient. 
2.   Inversion 

A  group  of  links  connected  together  with  no  links  rigidly  fixed 
is  called  a  kinematic  chain.   If  a  chain  is  arranged  so  that  the  motion 
of  one  link  causes  a  definite  motion  of  every  other  link  relative  to 
the  first  link,  it  is  called  a  constrained  chain.   Motion  of  one  link 
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of  the  chain  constrains  the  remaining  links  to  execute  determinate 
relative  motions.   A  constrained  chain  does  not  produce  any  definite 
absolute  motion.  When  one  of  the  links  is  fixed,  the  relative  motions 
of  the  chain  with  respect  to  that  link  become  absolute.   The  chain 
becomes  a  mechanism.   A  constrained  kinematic  chain  can  be  made  into 
a  mechanism  in  as  many  ways  as  it  has  links.   Some  of  these  mechanisms 
may  be  equivalent  due  to  conceptual  symmetry  of  the  chain. 

The  significance  of  inversion  is  that  the  motions  of  all  mem- 
bers of  a  chain  with  respect  to  each  other,  that  is,  their  relative 
motions  are  independent  of  which  member  serves  as  the  frame  or  fixed 
link.   It  will  be  shown  that  inversion  is  a  useful  tool  for  synthesis  of 
multi-link  mechanisms  because  it  allows  synthesis  techniques  for  simple 
mechanisms  (four-bar)  to  be  extended  to  the  design  of  more  complex 
mechanisms.   For  example,  each  six-link  chain  has  one  inversion  that  can 
be  synthesized  with  less  difficulty  than  the  others.   These  are  the 
Stephenson  3,  Fig.  (1-la)  and  the  Watt  2,  Fig.  (l-2a) .   Using  inversion, 
specifications  may  be  chosen,  e.g.,  for  a  Stephenson  2  mechanism  and  con- 
verted to  equivalent  synthesis  specifications  for  a  Stephenson  3. 
Synthesis  may  then  be  carried  out  and  the  resulting  Stephenson  3  mechan- 
ism reinverted  to  obtain  the  Stephenson  2  which  satisfies  the  chosen 
specifications. 

Although  the  equations  of  inversion  involve  only  angular  prop- 
erties of  links,  it  will  be  shown  in  Chapters  III  and  IV  how  inversion 
may  be  used  to  synthesize  mechanisms  for  specified  path  and  angular  prop- 
erties. A  Stephenson  3  mechanism  is  shown  in  Fig.  (l-3a)  for  two  positions 
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Figure  1-3  Inversion  of  a  Stephenson  Six-Link  Mechanism 
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of  the  chain.   The  exact  same  chain  is  shown  for  the  same  two  positions 
as  a  Stephenson  2  inversion  in  Fig.  (l-3b) .   The  angular  positions  and 
higher  order  time  derivatives  of  links  in  these  two  inversions  are 
related  by 

nj    nj    5  j    50 


i  .  =  6   -  9  n  =  link  number 

nj    nj    5j 


+  .  =  8  . 
nj    nj 


nj    rnj    T3j    r30 


j  =  position  counter  (1-1) 

n  =  initial  position  of 

link  5  in  Stephenson  2 


nj    nj    3j  ifiof,  ■  initial  position  of 

link  3  in  the 


Stephenson  3 


"J    nj    3j 


Note  that  9,  appears  because  link  5  is  fixed  in  the  Stephenson  3 
and  t(>3  appears  because  link  3  is  fixed  in  the  Stephenson  2.   Alterna- 
tive choices  for  fixed  links  are  2  and  4.   In  this  case  the  equations 
would  contain  8   and  t|i,  .   A  general  expression  for  any  chain  can  be 


n  =  link  number 

j  =  position  counter 

k  =  derivative  counter 

r  =  link  which  is  fixed  in     (1-2) 

mechanism  with  angles 

*  . 
nj 

s  =  link  which  is  fixed  in 

mechanism  with  angles 


ten 

(k) 

(k) 

(k)           (k) 

=      9    .    ■ 

-    e  .  +    8 

rj             rO 

(k) 

(k) 

(k)          (k) 

e  .  = 

=     *    .   ■ 

'sj          vs0 
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When  n  =  r  or  n  =  s,  Eq.  (1-2)  be 


and 


>l>    .   =  8  .  =  constant 


8„ .  =  i|>  „  =  constant 
sj    sO 


which  shows  that  r  is  the  fixed  link  in  the  mechanism  with  angles  ii 
and  s  is  the  fixed  link  in  the  mechanism  with  angles  9  . .   All  angles 
are  measured  from  the  same  initial  absolute  reference. 

The  computer  program  MECSYN  (Chapter  V)  determines  the 
inverted  specifications  for  a  general  three-link  chain  using  the  above 
equations  and  an  analysis  procedure  to  determine  path  properties.   This 
tool  is  used  to  synthesize  a  large  class  of  multi-link  mechanisms  by 
inversion  which  cannot  be  synthesized  by  other  closed  form  methods. 
Details  of  this  procedure  may  be  found  in  Chapters  III  and  IV  or 
References  [14,15,19]. 

3.   Path  Cognates 

Planar  four-bar  mechanisms  exhibit  a  property  known  as  the 
Roberts-Chebyshev  theorem: 

Three  different  planar  four-bar  mechanisms  will  trace  identical 
coupler  curves.   Mechanisms  related  by  this  property  are  known  as 
cognate  mechanisms,  Fig.  (1-4).   Samuel  Roberts  presented  a  proof  of 
this  theorem  in  1875,  while  Chebyshev,  using  an  entirely  different 
approach,  presented  his  proof  in  1878. 

Although  this  property  is  well  known,  some  extremely  important 
interrelations  within  this  set  of  mechanisms  are  not  widely  appreciated. 
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It  is  the  property  of  the  assembly  that  coupler  links  3,  6,  9  are 
directly  similar,  that  four-link  loops  I,  II,  III  are  parallelograms 
and  that  the  angular  displacement,  velocity,  acceleration,  and  all 
higher  order  properties  of  link  sets  {9,2,7},  {10,3,5}  and  {8,4,6} 
are  identical. 

Due  to  these  properties  of  cognate  mechanisms,  a  four-bar  may 
be  synthesized  for  path  and  angular  specifications  of  a  coupler  link 
(3  in  Fig.  (1-4))  giving  two  additional  four-bars  (8,9,10,1  and  5,6,7,1) 
whose  input  links  (10,5)  satisfy  the  desired  angular  specifications, 
while  the  common  coupler  point  C  satisfies  the  desired  path  specifica- 
tions.  The  result  is  a  known  coordination  of  a  crank  with  the  motion  of 
a  coupler  path  point.   This  tool  has  been  extended  to  aid  in  the  solu- 
tion of  many  multi-link  synthesis  problems,  especially  the  coordination 
of  a  crank  with  the  complete  motion  of  a  body  (path  and  angular  specs) 
[14]. 

The  analytical  formulation  for  the  calculation  of  four-bar  path 
cognates  may  be  found  under  the  heading  subroutine  COG  in  Chapter  V. 

4.   Angular  Cognates 

This  property  has  been  described  (by  Hain)  as: 

The  links  of  a  four-bar  mechanism  may  be  rearranged  in 
any  order  without  affecting  their  angular  motion. 

Note  that,  in  general,  the  same  coupler  point  location  for  two 

angular  cognate  four-bars  will  not  have  the  same  path  properties  as 

with  path  cognates. 
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Point  A  on  dyad  (1,2),  Fig.  (1-5),  traces  a  path  while  link  2 
undergoes  angular  motion  6.   A  parallelogram  (1,2,3,4)  is  formed  hy 
adding  dyad  (3,4)  to  dyad  (1,2).   Since  this  is  a  parallelogram,  link  4 
always  remains  parallel  to  link  2  for  all  motion  of  point  A  and  the 
angular  motion  <j>  of  link  4  is  related  to  link  2  by  a  constant  phase  angle. 
*■»  9  -  6   +  180° 


(k)   (k) 

*  =  B      for  all  k. 

Since  links  2,4  and  links  1,3  are  equal  in  length,  this  is  equivalent 
to  interchanging  links  1  and  2.   Applying  this  to  four-bar  (1,2,3,4) 
in  Fig.  (1-4),  first  to  dyad  BAO  and  then  to  dyad  ABO,  demonstrates 
the  property  of  angular  cognates  for  a  four-bar.   However,  this  property 
exists  for  a  single  dyad  alone  as  shown  above.   The  importance  of  this 
property  is  that  any  mechanisms  composed  of  dyad  and  four-bar  chains  can 
be  synthesized  for  angular  specifications  of  coupler  planes  and  have 
those  angular  properties  transformed  to  an  input  or  output  link  using 
angular  or  path  cognates,  respectively.   This  capability  is  essential 
for  the  synthesis  of  mechanisms  given  dynamic  specifications  as  will  be 
shown  in  Chapter  II. 
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Figure  1-5  Angular  Cognates 


CHAPTER  II 
THEORY  OF  DYNAMIC  SPECIFICATION  SYNTHESIS 

A.   Introduction 

The  coordination  of  geometric  motion  states  has  been  of  interest 
to  researchers  in  mechanisms  ever  since  the  time  of  Releaux.   In  most 
cases  the  motion  is  considered  to  be  dependent  on  a  single  position  param- 
eter; i.e.,  the  rotation  angle  8  for  an  input  crank.   The  output  state 
may  be  considered  as  being  functionally  related  <j>  =  f(e)  to  the  input 
parameter  by  means  of  higher  order  properties  <j>' ,<j>", . . .  of  this  function 
for  a  given  value  of  9.   This  type  of  motion  is  widely  considered  as  the 
specification  of  infinitesimally  separated  positions  (ISP)  for  the  geo- 
metric motion.   If,  however,  both  input  (q)   and  output  (((;)  motions  take 
on  finite  increments  A9  and  A$»  then  the  specification  is  considered  as 
giving  finitely  separated  positions  (FSP)  for  a  sequence  of  positions 
I   =  0,1,2,...  .   Note  that  for  FSP,  neither  parameter  9,$  need  be  consid- 
ered as  independent  to  specify  the  motion. 

To  illustrate  a  motion  involving  two  functional  constraints, 

consider  the  motion  (U,V)  of  the  point  A(u,v)  in  the  moving  plane  E  where 

9  is  the  independent  parameter  (see  Fig.  2-la) .   Given  the  functional 

dependence  x  =  f(0),  y  =  f (6) ,  due  to  the  mechanism,  it  is  obvious  that 

U  -   u  cos  9  -  v  sin  6  +  x(9) 

(2-1) 
V  =  u  sin  e  +  v  cos  9  +  y(e) 
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Figure  2-1   Fixed  and  Moving  Coordinate  Systems 
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which  is  the  explicit  coordination  of  the  motion  of  point  A  with  the 
parameter  9.   For  ISP,  one  must  have  knowledge  of  the  higher  order 
properties  x' ,x", . . .y ' ,y". . .  with  respect  to  the  independent  parameter 
0  to  obtain  U' ,U", . . .V , V", . . .  by  differentiating  Eq.  (2-1)  with  respect 
to  9.   The  derivatives  of  x(6),y(9)  are  the  constraint  conditions  which 
must  be  provided  by  the  mechanism. 

Now  suppose  the  motion  U,V  of  point  A  is  to  be  coordinated  with 
a  crank  parameter  <J>.   This  motion  involves  three  functional  constraints 
x  =  f($),  y  =  f(tf>)  and  6  =  f(<f>).   For  ISP,  it  is  necessary  that  the 
motion  be  specified  in  terms  of  higher  order  properties  of  these  func- 
tions with  respect  to  cf>.   For  the  system  being  driven  by  the  crank 
parameter  <f>,  an  important  dynamic  coordination  problem  evolves  when 
<J>  experiences  a  general  input  $  =   f(t)  with  a  given  set  of  higher  order 
properties  $,$,...  . 

Dynamic  state  of  a  system  is  usually  defined  as  the  smallest 

collection  of  numbers  which  must  be  specified  at  time  t  =  t  in  order  to 

r  o 

be  able  to  predict  uniquely  the  behavior  of  the  system  for  any  time 

tit  for  any  input  belonging  to  the  given  input  set  provided  that 

every  element  of  the  input  set  is  known  for  t  >  t  . 

o 

In  the  synthesis  method  considered  here,  in  some  cases  of  the 
motion  specification,  more  than  the  dynamic  state  is  actually  prescribed. 
Because  of  this,  the  term  dynamic  specification  or  dynamic  specs  is  used 
to  indicate  the  prescribed  set  of  conditions  for  synthesis. 

For  example,  if  the  inertia  force  and  inertia  torque  are  prescribed 
for  an  object  at  a  given  position  having  mass  m  and  rotational  moment 
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of  inertia  I  (the  object  is  moving  in  coplanar  motion),  the  linear  and 
angular  velocities  and  acceleration  of  the  mass  center  may  be  computed. 
(Note  that  the  object  need  not  be  undergoing  coplanar  motion  but  may  have 
general  spatial  motion.   The  dynamic  specs  synthesis  is  still  valid  but 
must  then  be  modified  to  use  the  constraint  function  and  subsequent 
solution  algorithms  to  reflect  the  type  of  spatial  constraint  pairs 
used.)   Higher  order  properties  which  influence  motion  constrained  sys- 
tems may  also  be  prescribed,  such  as  jerk  (time  derivative  of  acceler- 
ation) and  jounce  (time  derivative  of  jerk).   Synthesizing  a  mechanism 
which  satisfies  these  specifications  while  driven  by  a  prescribed  input 
determines  the  dynamics  of  the  mechanism  at  that  position.   The  pre- 
scribed inertia  force  and  torque  are  achieved. 

The  solution  of  this  problem  means  that  the  motion  of  the  moving 
plane  E(u,v)  as  x  =  f(t),  y  =  f(t),  8  =  f(t)  can  be  coordinated  with 
a  general  input.   The  objective  here  is  to  give  a  fundamental  procedure 
for  treating  the  synthesis  of  analytically  coupled  systems  for  the  4-, 
5-,  and  6-link  planar  mechanisms. 
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B.   Previous  Synthesis  Philosophy 

If  the  motion  depicted  in  Figure  (2-la)  is  considered  as  con- 
strained, then  6  may  be  considered  as  the  independent  motion  parameter 
and  x  =  f(6),  y  =  f(6)  become  the  conceptual  representation  of  the  phys- 
ical constraint.   The  fundamental  form  of  the  physical  constraint  is 
the  crank  where  the  pin  joint  A(u,v)  in  E  moves  on  a  circular  arc  in 
the  reference  plane  E.   It  is  possible  to  consider  a  general  functional 
constraint  F(U,V)  =  0  for  the  motion  of  point  A  which  has  the  internal 
coefficients  0^.   Using  the  transform  given  in  Eq.  (2-1),  the  conceptual 
constraint  becomes 

F  =  f(u,  v,  x,  y,  6,  QQ,  0^,  Q2>...)  =  0  .  (2-2) 

For  FSP  motion  specifications  {x.,y.,0.}  for  the  motion  of  plane  E,  the 
problem  becomes  defined  by 

Fj  =  f(u,  v,  Xy    y.,    By    QQ,  Q1>  Q2,...)  =  0  (2-3) 

J  -  0,1,2,... 
where  the  unknowns  are  the  coordinates  (u,v)  and  coefficients  Q.. 
For  ISP,  the  geometrical  conditions  on  the  motion  are 

de 

k  =  0,1,2,... 
for  a  given  value  of  9  where  the  unknowns  are  the  same  as  for  FSP.    s. 

The  concept  of  MSP  implies  the  mixing  of  these  conditions  to 
provide  a  means  of  more  generally  specifying  the  motion.   If  the  symbol 
PP  denotes  2  ISP,  and  P-P  denotes  2  FSP,  then  the  following  combinations 


Fk  ' k  tf(u'  v>  x>  y>  9>  V  V  Q2»"->}  "  °        (2-4) 
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PPPPP,  PPPP-P,  PPP-PP,  PPP-P-P 
PP-PP-P,  PP-P-P-P,  P-P-p-p-p 
represent  the  potential  for  7  distinct  motion  problems.   The  goal  should 
be  to  establish  a  uniform  body  of  analytics  which  treats  all  of  these 
problems  simultaneously;  i.e.,  the  basis  for  the  MSP  concept.   If  I 
counts  all  MSP,  j  counts  all  FSP  and  k  counts  all  sequential  ISP  for  a 
distinct  position  j ,  then  the  case  symbol  is  sufficient  t.o  denote  the 
kind  of  motion  desired.   For  example,  for  case  PPP-PP,  the  following 
counters  apply 

I  0     1     2     3     4 

j     0     0     0     1     1 

k     0     1     2     0     1    . 

Generally,  the  counter  I   is  understood  to  represent  the  other  counters 
whenever  the  case  symbol  is  given. 

C.   Statement  of  Dynamic  Synthesis 

The  majority  of  previous  synthesis  problems  have  been  referred 
to  as  kinematic  synthesis  solutions  when  in  fact  they  have  not  included 
specifications  for  kinematics  in  the  strictest  definition  of  the  term 
(such  as  velocity,  acceleration,  jerk,  and  jounce).   The  objective  has 
been  to  synthesize  mechanisms  to  satisfy  purely  geometric  kinematic 
properties  (such  as  the  tangent,  curvature,  inflexion,  etc.,  of  a  given 
function).   Here  the  goal  is  to  synthesize  mechanisms  for  prescribed 
time  dependent  properties  in  the  form  of  velocity,  acceleration,  etc., 
which  along  with  specified  mass  and  rotational  inertia  make  up  the 
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inertial  dynamics  of  the  coupler  link.   The  problem  will  be  formulated 
in  terms  of  six-link  mechanisms,  Fig.  (2-2),  for  two  general  types  of 
solvable  problems.   Those  in  Fig.  (2-3)  involve  the  central  concept  of 
path  generation  with  auxiliary  angular  motions  being  coordinated. 
Those  in  Fig.  (2-4)  deal  completely  with  the  coordination  of  angular 
motions.   These  solvable  problems  have  been  treated  in  complete  detail  in 
Ref.  [15].   The  level  of  the  motion  specification  depends  on  the  problem 
but  in  general  involves  partial  dynamic  specification  (the  dynamics  related 
to  a  point  or  the  angular  dynamics)  of  several  links  in  addition  to  the 
input.   The  prescribed  quantities  for  these  two  types  of  problems  are: 
Type  I   (See  Fig.  (2-3)) 


Prescribed  Dynamic 

Specifications  (Add  I,M)  of    (2-5) 

a  Moving  Body 


ar 

x.,i.,...x. 

(k) 

> 

yj-V-yj 

(k) 

e.,e e. 

(k)  1 

♦j.tj.  —  tjj 

Type  II      (See  Fi 

I-    (2-4)) 

(kf 

Vsr---ej 

) 

(k) 

iK.lJ.f.tJl. 
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Figure  2-3  Solvable  Problems  -  I 
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where  the  counters  (k,j)  for  position  2  depend  on  the  combination  (case) 

,  00  dkx 

of  the  multiply  separated  positions  and  x  =  — r-  . 

dt 

To  satisfy  these  specifications  for  a  six-link  mechanism,  it  is 

necessary  to  find  two  or  more  points  in  the  appropriate  coupler  links 

that  satisfy  the  following  conditions: 

1.  The  points  must  travel  on  a  circular  arc  in  the  reference 
system  E(U,V). 

2.  The  points  experience  a  time  dependent  motion  which 
assures  that  the  dynamic  specs  of  the  link  are  satisfied. 
After  synthesis,  this  condition  is  met  by  the  coordination 
of  the  mechanism  with  the  prescribed  input. 

Using  the  same  philosophy  of  Section  B  for  kinematic  synthesis,  the 

dynamic  synthesis  constraint  function  for  circular  motion  of  a  point  can 

be  used  to  establish  the  conditions  to  be  met;  i.e., 

{G(x(t),    y(t),    6(t),    *(t)))|  =   0 

|t=t£ 

^{G(x(t),    y(t),    B(t),    <Kt))}|  =  o 


'U 


(2-7) 


dk 


-iMG(x(t),   y(t),    6(t),    d)(t))}|  =  0 

dtk  lt-tt 


where  the  symbols  x,  y,  0,  0,  may  correspond  to  those  in  Fig.  (2-lb). 
Satisfaction  of  these  conditions  for  mixed  FSP  and  ISP  will  be  called 
dynamic  specs  multiply  separated  positions  (DMSP) . 

The  author  has  shown  in  Refs.  [14,151  that  this  constraint  link 
formulation  or  curvature  relation  can  be  decoupled  for  purposes  of  design 
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by  means  of  the  path  point  for  which  the  dynamic  specifications 
are  made.   This  process  allows  a  modular  design  of  five-link,  six-link, 
and  more  complex  mechanisms  if  desired.   Thus,  the  four  tools  of  syn- 
thesis are  applied  to  at  most  a  four-link  chain  or  frequently  only 
a  dyad.   Two  examples  will  be  given  here  which  illustrate  this  modular 
approach  for  the  Type  I  and  Type  II  dynamic  specs  synthesis. 

Conceptually,  it  is  possible  to  decouple  the  design  procedure 
into  two  or  more  modules  which  allows  the  formulation  of  two  distinct 
sets  of  constraint  conditions  in  the  form: 

,k 
-S-  {G(x(t),  y(t),  6(t))}|      =0  (2-8) 


dt  I  ft 


I 


and 


,k 
-2-  {G(x(t),  y(t),  4>(t))}|     =0.  (2-9) 

dtk  |«f 

Example  1  is  associated  with  the  Stephenson  3  for  the  solvable  Type  I 

problems  depicted  in  Fig.  (2-3).   (See  Refs.  14, 15, 19].)   The  problem 

specification  for  the  problem  shown  in  Fig.  (2-5)  can  be  established  by 

the  following  steps: 

Step  (1)   Choose  the  dynamic  specs  for  the  coupler  link  by  the  specs 
(k)    (k)  (k) 

{  x.(t),  y.(t)}  at  point  C  and  {  8.(t)}  for  the  angular 

motion  of  that  link.   Find  one  constraint  link  (link  1)  by 

using  the  dynamic  curvature  transform  respresented  by 

Eq.  (2-8).   (See  Fig.  (2-5a).) 


The  path  point  could  be  a  pin  joint  (such  as  C)  or  some  other 
associated  point  in  the  system. 
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Step  (2)   Use  the  same  specs  for  point  C  along  with  the  desired 
(k) 
input   set   $.(t)   for  a  coupler  link.   Find  two  con- 
straints (links  2  and  3)  by  the  curvature  transform  related 
to  Eq.  (2-9)  to  form  a  four-bar  linkage.   Taking  the  path 
cognates  of  this  four-bar  results  in  a  four-bar  with 

00 

cranks  2'  and  3',  one  of  which  contains  the  set   *.(t) 

as  the  input.   (See  Fig.  (2-5b).) 

Step  (3)   Steps  (1)  and  (2)  have  created  systems  which  have   common 

00     (k) 
specs  at  point  C  {  x.(t),   y.(t)}.   Join  these  at  point  C 

to  obtain  the  desired  Stephenson  3  six-link  mechanism 

(Fig.  (2-5c).) 

Example  2  uses  the  Stephenson  3  to  illustrate  the  type  of  synthesis 

process  involved  with  the  Type  II  problem  (see  Ref.  [15])  and  depicted 

in  Fig.  (2-4).   It  involves  the  following  steps: 

(k)  (k) 

Step  (1)   Choose  the  angular  dynamic  soecs  {  6 .  (t) ,  Oi  (t) }  for  two 

J  J 

adjacent  links  of  chosen  length  H    ,    I     and  location 

A(U  ,V  )  in  the  reference  system  £(U,V)  (see  Fig.  (2-6a)). 

Analyze  this  dyad  to  determine  the  dynamic  specs  of  point  C 


(k)      dk 

x.(t)  =  — -  {£  cos  9  +  I      cos  ii   +  U  }| 

3  dtk   1  2  a  t-t 

'    J 

(k)      dk 

y.(t)  =  — r-  U  sin  6  +  £-  sin  ill  +  V  }  | 
J      dtk   l  2  a  |  t-t 

(k) 
Step  (2)   Specify  the  input  set   <|> .  (t)  desired  to  drive  the  system 

along  with  the  specs  of  point  C  to  create  the  dynamic  specs 
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-oC(Xj(t),y/t)) 


C(Xj(t),yj(t)) 


b. 


f  cAm&ti) 


Figure  2-5  Type  I  Synthesis  Problem 
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a. 


CttM.jM) 


Type  2-6  Type  II  Synthesis  Problem 
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for  a  coupler  link.   Using  the  curvature  transform  repre- 

00    00    GO 

sented  by  Eq.  (2-9)  for  the  specs  {  x.(t),   y.(t),   <t.(t)}, 
find  two  constraint  links  to  form  a  four-bar  linkage. 
Taking  the  path  cognate  of  this  linkage  results  in  a  four- 
bar  with  cranks  3'  and  4',  one  of  which  contains  the  set 
(k) 
<p .  Ct )  as  shown  in  Fig.  (2-6b)  . 

(k) 
Step  (3)   Steps  1  and  2  contain  point  C  with   common  specs  {  x  (t) , 
(k)  J 

y. (t)}.   Joining  these  systems  as  shown  in  Fig.  (2-6c) 

gives  the  desired  Stephenson  3  mechanism. 
In  this  fashion  more  than  35  distinct  six-link  synthesis  problems  accord- 
ing to  the  step-by-step  procedures  of  Refs.  [14,15,19]  are  decoupled  into 
component  problems  of  low  complexity.   The  examples  here  illustrate  the 
least  difficult  of  these  procedures.   Chapter  IV  will  consider  more 
involved  procedures  to  completely  validate  the  generality  of  this  section. 
This  generality  is  distinguished  from  the  author's  previous  work  by  the 
dynamic  curvature  transform  which  will  be  presented  in  the  next  section. 
All  the  other  tools  of  synthesis  (inversion,  angular  cognates,  and  path 
cognates)  maintain  their  previous  significance. 

D.   Formulation  of  the  Dynamic 
Curvature  Transform 

The  constraint  conditions  used  to  synthesize  mechanisms  depends 
on  a  circular  constraint  which  can  be  represented  as 

g(U,V)  =  QQ(U2+V2)  +  2Q1U  +  2Q2V  +  Q3  =  0  (2-10) 
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which  for  dynamic  constraint  conditions  requires  that 

dk 
Sfc k  (Q0(U2+V2)  +  2Q  U  +  2Q  V  +  Q  }|     =0  (2-11) 

dt  |t=t* 

for   I   =  0,1,2,...  . 
Subtracting  the  condition  for  I   =  0  gives  the  generalized  circle  con- 
straint function 


Gl   =  Vg0  =  ~ k  {Qo(u   +V  }  +  2Q1U  +  2Q2V}I     =  °         <2"12> 
dt  I t=t. 

for   l   =  1,2,3,... 

which  removes  the  coefficient  Q3  from  further  consideration.   Using  the 
coordinate  transformation  (U,V)  -*  (u,v)  given  in  the  introduction  trans- 
forms this  system  into  the  form 

Gz  =  DcnQo  +  DnuQo  +  D2£v%  +  d3ji(uQi4^2) 

+  D4t(uQ2-vQ1)  +  D5iQl  +  D6iQ2  =  0  X2-13) 

I   =  1,2,3,4 

where  I   is  the  counter  corresponding  to  the  (k,j)  indices  for  the  case 

of  DMSP.   The  fundamental  coefficients  D  .  may  be  considered  as  the 

dynamic  circular  constraint  coefficients  and  are  listed  for  all  possible 

(k,j)  values  in  Table  2-1.   To  establish  the  constraint  condition  G  , 

one  chooses  the  U^   from  the  table  according  to  the  indices  (k,j)  for 

the  value  of  the  counter  £.   These  coefficients  are  completely  general 

and  are  applicable  for  any  location  of  the  reference  systems  E(u,v)  and 

£(U,V).   All  the  D^  have  been  tabulated  in  a  single  computer  program. 

Since  all  cases  of  DMSP  are  represented  simultaneously  by  the  D 

m£' ' 

a  considerable  savings  in  programming  complexity  results.   In  addition, 
all  synthesis  equations  can  be  formulated  explicitly  in  terms  of  the 
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D   making  it  possible  to  create  a  single  body  of  analytics  to  cover  all 
m£ 

cases.   This  means  that  one  (not  seven)  computer  program  is  sufficient 
to  treat  DMSP.   The  modular  synthesis  approach  based  on  sub-unit  link 
groups  (the  dyad  and  the  four-bar)  allows  synthesis  of  many  different 
types  of  planar  mechanisms  in  terms  of  a  relatively  small  computer  pro- 
gram rather  than  the  cumbersome  alternative  of  a  program  for  each 
mechanism. 

Equation  (2-13)  consists  of  four  equations  which  are  nonlinear 
in  the  circle  point  coordinates  (u,v)  and  the  centerpoint  coordinates 

Ql         Q2 
(U  =  -  —  ,  V  =  -  -— )  of  the  circle  on  which  the  circle  point  moves 

(see  Eq.  (2-10)).   To  obtain  a  solution,  a  mapping  into  a  new  set  of 

variables  Z   (m  =  0,1,... 6)  is  defined  by 
m 

ZQ  =  QQ        Z3  =  UQ1  +  VQ2        Z5  =  Ql 

Zl  =  uQ0       Z4  =  uQ2  "  vQl        Z6  =  Q2  (2_U) 

Z2  =  ^0 
which  means  that  Eq.  (2-13)  can  be  written  as  a  linear  system  of  equa- 
tions in  the  form 

6 

I     D  .  Z  =  0  I   =  1,2,3,4  .  (2-15) 

„  mJ.  m  ' 

m=0 

Two  additional  unknowns  have  now  been  introduced  to  allow  the 
formulations 


Z0Z3  =  Z1Z5  +  Z2Z6 


Z0Z4  =  Z1Z6  "  Z2Z5 


(2-16) 


Equations  (2-15)  and  (2-16)  constitute  six  equations  in  the  six  var- 
iables Z _/Z  ,  m  =  1,2,... 6  yielding  four  sets  of  desired  solutions 
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Center 
Z  ] /"Points 


Dynamic 
Specifications 
Constraint      (2-17) 
Links 


These  pairs  of  points  could  be  considered  to  be  generalized  Burmester 
pairs  for  dynamic  synthesis. 

The  actual  solution  of  Eqs.  (2-15)  and  (2-16)  is  established  by 
first  solving  the  linear  system  of  Eq.  (2-15)  for  the  unknowns  Z  ,  Z  , 
Z5'  Z6  in  terms  of  the  unknowns  Z  ,  Z  ,  Z  in  the  form 


D31  D41  D51  D61 


D32  D42  D52  D62 


D33  D43  D53  D63 


D,.   D,  .   Dc,   D, , 
34   44   54   64 


Z3 
Z4 
Z5 
Z6 


D01  °11  D21 


D02  D12  D22 


D03  D13  D23 


V  D14  °24 


(2-18) 


B>Z  =  b  (2-19) 

where  the  j,k  correspond  to  the  l   counter  for  the  case  being  considered. 
When  ]P    is  non-singular,  the  solution  is 


and 


Z  =  Z  A  +  Z  B  +  Z.C    m  =  3,4,5,6 
m    0  m    1  m    2  m  ' 


|D|    l!l   D°£   A"ni'  Bm  =   "    |p  |       J.J.  D"   Am<l    ' 


(2-20) 


c„  -       -*-     Z     D„„    A 


d   „ 
m£  "        mJl 


28.     mJl 


JP   -     lldmt   ll>    AmP    is    the   cofactor  of   element   d      . 

mJl 
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WhenjJ)  is  singular  as  it  is  in  the  case  of  a  specified  motion  for 

a  point  along  a  line,  then  the  unknowns  Z,,  Z„,  Z„,  Z  are  eliminated 

12        3        4 

in  favor  of  ZQ,    Z  ,    Z    .      This  alternative  process    treats  most  special 


Eliminating  Zy    Z^,   Zy   Z&  between  Eq.    (2-16)    and  Eq.    (2-20)    gives 
two  general  conies   in   the   coordinates  u  =   Z   /Z     and  v  =  Z   /Z    : 

av2  +  buv  +  cv  +  du2+eu  +  f     =0 

(2-21) 
a'v2  +  b'uv  +  c'v  +  d'u2  +  e'u  +  f '  =  0 

a  =  C4,   b  -  C3+B4,   c  =  A4-Cl>   d  =  B3,   e  -  Aj-B^  f  =  -^ 

a'  =  Cj,   b'  =  B3-C4>   c'  =  A3-C2,   d'  =  -B4>  e'  =  B^,  V    =  A.,. 

Using  the  Bezout  Eliminant  [23]  for  Eq.  (2-21),  a  quartic  polynomial 
in  u  results  as 

V"  +  V3  +  E2u2  +  Exu  +  EQ  =  0  .  (2-22) 

If  (ab')  =  ab'  -  a'b,  the  coefficients  may  be  written  in  the  form 
E4  =  (ad')2  -  (ab')(bd') 

E3  =  2(ad')(ae')  -  (ab')[(be*)  +  (cd')]  -  (ac')(bd') 
E„ 

-  (ac')[(be')  +  (cd')] 
Vx   =  2(ae')(af)  -  (ac'H(bf')  +  (ce')]  -  (ab')(cf) 
E0  =  (af')2  "  Cac')(cf').  (2_23) 

The  common  root  v  is  given  by 

„  -[(ad')u2  +  (ae')u  +  (af')l 

[(ab')u  +  (ac')]       •  (2-24) 


-2  =  (ae')2  +  a(ad')(af)  -  (ab')t(bf')  +  (ce')] 
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The  Burmester  circle  points  for  the  dynamic  specs  problem  are  (u  ,v  ) 
corresponding  to  the  max  (n)  =  0,2,  or  4  solutions  of  Eq.  (2-22). 
The  associated  center  points  to  complete  the  data  in  E  for  the  con- 
straint pair  are 

U  -  -  (A.  +  Bcu  +  C,v  )  (2-25) 

n       5    5  n    5  n 

V  =  (A,  +  B.u  +  C,v  )    max  (n)  =  0,2,4  . 
n     o    b  n    b  n 

For  motion  of  a  point  along  a  line,  Eq.  (2-22)  reduces  to  a  cubic  and 
max  (n)  =  1  or  3. 

E.   Problem  Concepts 

Ordinary  geometric  or  kinematic  synthesis  may  be  considered  as 
a  special  case  of  this  formulation  and  can  easily  be  obtained  by  replac- 
ing t  by  8  (t  —  0)  and  by  requiring  that 

'-■  s-s-*.  $-■$-■    >>* 

at    do 
such  that  9  becomes  the  independent  parameter  in  the  constraint  function 

dk 

G£  =   k  {G(x(t)'  yW>  e(t))}|     =  0  (2-27) 

dt  I  t-h 

which  now  becomes 

dk 

G,  =— r{G(x(8),  y(0),  8)}|     =0   .  (2-28) 

da1"  I  e-e^ 

The  dynamic  coefficients  D  „  of  Table  2-1  automatically  reduce  to  the 
AfflJl  of  kinematic  synthesis  for  MSP  found  in  Ref.  [  9  ]. 

Another  useful  property  of  this  formulation  is  easily  obtained 
merely  by  changing  the  input  specifications  to  the  program.   This 
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involves  the  geometric  synthesis  for  MSP  of  complex  six-link  mechanisms. 
The  previously  mentioned  Type  I  and  Type  II  problems  now  become  speci- 
fied as 


Prescribed  Geometric  Specs 
of  a  Moving  Body 


Type   I 

¥«•  (f 

V-©, 

',"'■  (% 

\~$\ 

V«'  (f 

)r<\ 

Type   II 

•,<♦>•  (S 

v-®,} 

v-ev-($) 


Prescribed  Geometric 
Angular  Specs 


where  the  independent  parameter  $ .  takes  on  increments  for  FSP  only 
(i.e.,  <j>  =<(>,).   The  Type  I  problem  is  the  coordination  of  a  coupler 
path  x(<(j) ,  y(<|>)  with  a  coupler  angle  8(<(>)  (or  output  angle)  and  an 
input  angle  $.   The  Type  II  problem  is  the  coordination  of  three  angu- 
lar quantities  with  <f  the  input.   The  kinematic  problems  of  the  Type  I 
and  Type  II  specifications  which  were  not  previously  treated  by  the 
author  [15]  can  now  be  immediately  solved  in  the  following  manner. 

The  independent  parameter  t  may  be  considered  as  a  dummy 
parameter  and  replaced  by  an  input  angle  (t  *v   8)  while  reinterpreting 
Eq.  (2-7)  as 
,k 


{G(x(9),   y(9),  8,  (K9))  )| 


=  0 


do 


(2-29) 
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This  interpretation  is  used  in  modular  form  in  three  steps  similar  to 
the  use  of  Eqs.  (2-8)  and  (2-9)  in  Example  1  of  Section  C  to  design 
complex  six- link  mechanisms.   The  condition  following  from  Eq.  (2-8)  for 
the  first  step  becomes 

Hk 
-S-  {G(x(B),  y(8),  9)}      =0  (2-30) 

deK  I 9=ea 

which  reduces  to  an  ordinary  kinematic  (or  geometric)  synthesis  problem 
and  can  be  solved  in  the  usual  manner  (i.e.,  by  using  the  MSP  formula- 
tion in  Ref.  [  9  ]  or  the  complex  vector  loop  equation  in  Ref.  [24])  of 
by  using  the  "dynamic  MSP "  synthesis  program  with  the  specifications 

dV(t)    dkx.(9)   dky.(t)    dky.(6) 

J + J ,  J + 1_  (2-31) 

dtk      d6k      dtk      dB 

„  de   de   ,   d2e   d2e   d4e 

in  the  computer  program  input.   The  second  step  cannot  be  solved  by  the 
usual  kinematic  synthesis  methods  and  must  be  done  by  the  "dynamic 
MSP"   formulation.   In  Step  2,  Eq.  (2-9)  becomes 

-^j-  {G(x(6),  y(6),  <K6))}|  (2-32) 

d9  I  9=6^ 

which  is  distinct  from  Eq.  (2-30).   Since  t  is  replaced  by  9,(t  %   9), 
the  specifications  are  modified  to 


dx.(t)        dkx.(8) 

.1           .           .1 

dky.(t)        dky.(9) 

dk<t>.(t)      d%.(e) 

k                     k 
dt                   d8K 

dtk            dek 

k       +            k 

dt                   d9K 

where  the  indices  j,k  depend  on  the  counter  i   and  the  case  being  con- 
sidered. 
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Part  II  presents  two  practical  applications  of  dynamic  specs 
synthesis  (DMSP) .   The  first  example  utilizes  the  full  capabilities  of 
the  method  in  designing  a  Stephenson  3  and  a  Stephenson  2  for  the  same 
dynamic   specifications.       The  second  example  of  geometric  MSP  for 
kinematic  synthesis  as  described  in  terms  of  Eqs.  (2-30)  and  (2-32)  uses 
the  dynamic  specs  formulation. 

The  philosophy  of  dynamic  specs  synthesis  shows  that  a  large 

range  of  unsolved  problems  can  now  be  solved  for  not  only  coordinated 

geometric  specs  but  dynamic  specs   as  well.   Because  of  the  generality 

of  the  D  „,  a  single  body  of  analytics,  uniform  in  its  format  is  now 
m*. 

made  available  in  terms  of  a  single  computational  tool.   As  stressed  in 
Refs.  [15,13  and  extended  in  Ref .  [14],  the  principal  task  of  the  designer 
is  the  determination  of  the  proper  set  of  specifications  for  a  given 
mechanism  constraint  to  satisfy  the  posed  problem.   For  geometric  syn- 
thesis it  has  been  possible  to  establish  all  properties  of  the  system 
(such  as  pole  locations,  generalized  inflection  circles,  cubic  circle 
point  curves,  Burmester  point  coordinates,  etc.)  in  terms  of  the  A . 
Because  the  concept  of  dynamic  synthesis  is  a  logical  extension  of 
geometric  synthesis,  it  is  clear  that  the  D   will  allow  the  same 
explicit  formulation. 

The  author  has  developed  a  complete  and  easily  used  computational 
tool  for  the  expanded  concept  of  DMSP.   This  computer  program  is  called 
MECSYN  and  will  soon  be  made  available  with  a  complete  user's  manual  to 
assist  the  designer  in  its  operation.   MECSYN  is  a  non-iterative,  effi- 
cient program  containing  the  four-  and  five-position  dynamic  specs  MSP 
curvature  transform  for  all  seven  cases  in  addition  to  path  cognates, 
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angular  cognates,  inversion,  link  length  ratio  and  Grashof  criteria 
calculation.   All  special  cases  are  included  such  as  straight  line 
motion.   MECSYN  also  contains  an  analysis  package  which  allows  the 
designer  to  investigate  the  usefulness  of  the  synthesized  mechanisms. 
Use  of  MECSYN  will  be  described  in  Chapter  V   along  with  a  listing 
of  the  current  version. 


CHAPTER  III 


DYNAMIC  SPECS  SYNTHESIS  OF 
FOUR-,  FIVE-  AND  SIX-LINK  MECHANISMS 


In  this  chapter,  the  synthesis  procedures  will  be  developed  and 
illustrated  by  examples  for  representative  mechanisms  (as  shown  in 
Fig.  (2-3)  and  Fig.  (2-4)  and  cases  utilizing  the  four  tools  of  dynamic 
synthesis  as  described  in  Chapters  I  and  II.   First,  a  four-bar  dynamic 
specs  "function  generator"  is  designed  for  Case  2  (PPPP-P)  using  both 
angular  and  path  cognates.   Next,  a  Stephenson-2  six-link  mechanism  is 
designed  to  coordinate  the  angular  dynamic  specs  of  three  links  for 
Case  1  (PPPPP).   Finally,  a  geared  five-link  mechanism  is  designed  to 
coordinate  the  input  with  the  dynamic  specs  of  a  coupler  link  for 
Case  4  (PPP-PP) .   The  examples  in  this  chapter  are  verified  by  two  dis- 
tinct analysis  methods,  one  numerical  and  one  closed  form  and  non- 
iterative.   The  synthesis  procedures  for  the  application  examples  in 
Chapter  IV  are  verified  by  a  third  distinct  non-iterative  analysis 
method.   Also,  the  examples  in  Chapter  IV  demonstrate  the  remaining  sig- 
nificant case  of  dynamic  specs  synthesis,  i.e.,  Case  3  (PPP-PP).   The 
other  cases  are:   Case  7  (P-P-P-P-P)  which  is  purely  geometric  kinematic 
synthesis  by  definition,  since  it  involves  specification  of  finitely 
separated  positions  (FSP)  only;  Case  5  (PP-PP-P)  and  Case  6  (PP-P-P-P) 
which  involve  specification  of  finitely  separated  positions  (P-P)  or 
velocity  distributions (PP) .   Since  geometric  kinematic  synthesis  is 
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a  special  case  of  dynamic  specs  synthesis  all  seven  of  these  cases  are 
significant  for  the  specification  of  geometric  kinematic  properties 
such  as  position,  tangent,  curvature  and  rate  of  change  of  curvature 
referred  to  as  multiply  separated  positions  (MSP). 

A.   Four-Bar  Dynamic  Specs  "Function  Generator" 

The  name  "function  generator"  is  the  classical  name  used  for 
a  mechanism  whose  output  link  position  parameter  approximates  a  pre- 
scribed function  with  respect  to  its  input  link  position  parameter, 
i.e.,  pure  geometric  kinematic  synthesis.   Regardless  of  what  the 
input  dynamics  are*  the  mechanism  still  generates  the  same  geometric 
function.   In  the  context  of  dynamic  specs  synthesis,  it  is  a  misnomer. 
In  fact,  a  four-bar  synthesized  to  coordinate  the  input  with  the 
dynamic  specs  of  the  output  approximates  two  functions,   output  vs. 
time  and  input  vs.  time.   This  again  illustrates  the  dependence  of  the 
dynamic  specs  synthesis  on  the  appropriate  time  dependent  input. 
The  resulting  mechanism  always  satisfies  the  same  output  vs.  input 
function  regardless  of  what  time  dependent  input  is  provided.   This 
follows  since  the  mechanism  function  is  only  dependent  on  the  invariant 
link  dimensions.   The  name  function  generator  is  used  here  to  draw  atten- 
tion to  the  contrast  between  classical  geometric  kinematic  synthesis 
and  dynamic  specs  synthesis. 


* 

Input  and  output  refer   to  links  connected  to  ground  as  opposed 
to  coupler  links — not  connected  directly  to  ground. 
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C*=( -1.695,-1.301) 
C»  (1.16, .837) 
0  •  M.535,-1.464) 
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CASE  2 

pos.  der. 

e        <t 

0      0 

45°               0° 

0      1 
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0      2 
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0      3 

15  rad/»«'  8  rt>d/j»c3 

1       0 

85°             20° 

Figure  3-1  Four-Bar  Function  Generator  by  Angular  Cognates 
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Rather  than  using  inversion  to  synthesize  the  four-bar  func- 
tion generator  as  has  been  done  in  the  past,  a  new  method  will  be  used 
as  presented  by  the  author  and  Tesar  (Ref.  [15])  using  angular  cognates. 
The  advantage  of  this  method  (which  is  less  difficult  conceptually  than 
inversion)  is  that  it  gives  an  alternate  set  of  solutions.   This  is  not 
the  case  with  the  inversion  method,  since  the  angular  specifications 
for  one  of  the  inverted  links  is  a  function  of  both  input  and  output 
specifications . 

The  case  treated  here  involves  the  specification  of  velocity, 

acceleration  and  jerk  at  the  first  position  and  a  second  position. 

Refer  to  Fig.  (3-1)  in  the  following  design  steps: 

Step  1.   Prescribe  the  angular  dynamic  specs  of  two  adjacent  links 
(k)     (k) 
(4,3'),  as  {  6.(t),   <j>.(t)}  (input  and  output)  and  also 

the  location  of  points  A,  B  and  C'  as  shown  in  Fig.  (3-1): 

(k)  (k) 


SL  j  k 

8 . (t)  -  input 

((».  (t)  -  output 

0  0  0 

45° 

0° 

10  1 

5  rad/sec 

2  rad/sec 

2  0  2 

2 
10  rad/sec 

4  rad/sec 

3  0  3 

3 

15  rad/sec 

3 
8  rad/sec 

4  10 

85° 

20° 

A  =  (0,  0),     B  =  (0,  1),     C'  =  (1.5,  1). 
Note  that  the  prescribed  initial  angle  6   of  the  input  will  be  the 
actual  resulting  initial  angle.   In  general,  this  is  not  true,  since 
usually  the  angular  differences  between  finitely  separated  positions 
are  preserved.   Using  MECSYN,  the  equivalent  synthesis  specifications 
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for  point  C'(x,y)  are  found  by  analyzing  dyad  ABC',  represented  by 

(k)       k  _        

X.(t)  =  — r  {AB  cos  e  +  BC'  cos  $  +   X  }, 

J      dtk  a  lt=t. 

J 

(3-1) 

(k)       k   _        

y  (t)  =  — -  {AB  sin  a  +  BC'  sin  *  +  Y  } 

j      dt  a  | ft. 

J 

Step  2.   Use  the  dynamic  specs  curvature  transform  of  MECSYN  as 

formulated  in  Chapter  II  with  the  resulting  dynamic  specs 

(k)      (10      (k) 
from  Step  1,  {  X.  (t) ,   y.(t),    <f>.(t)}  . 


1 

j 

k 

(k) 

x.(t) 

(k) 
y.(t) 

(k) 
*.(t) 

0 

0 

0 

1.5  in 

1  in. 

0° 

1 

0 

1 

-5  in./sec 

6  in/sec 

1  rad/sec 

2 

0 

2 

-37  in./sec 

-13  in./sec 

4  rad/sec 

3 

0 

3 

3 

-52  in/sec 

3 
-260  in/sec 

3 

8  rad/sec 

4 

1 

0 

0.5931  in. 

1.5798  in. 

20° 

One  of  the  solutions  given  by  MECSYN  is  the  originally  chosen 
constraint  link  AB. 

One  of  the  remaining  solutions  is  the  constraint  link  C"D  forming 
four-bar  ABC'D  shown  in  Fig.  (3-la).   Although  at  first  glance  this 
appears  to  be  a  crank-rocker,  MECSYN  indicates  it  is  a  double-rocker 
mechanism  with  link  length  ratio  15.54. 

C"  =  (-1.695,   -1.301),    D  -  (-1.535,  -1.464). 
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■SteP  3-  The   angular  cognate  of  dyad  DC"B  is  calculated  to  give 
the  angular  cognate  point,   C  =  (1.16,  .837). 
Point  C  completes  the  parallelogram  DC'BC  shown  in 
Fig.  (3- la).   Link  DC  now  has  the  desired  output  angular 
dynamic  specs.   Figure  (3-lb)  shows  the  two  prescribed 
finitely  separated  positions  of  the  resulting  four-link 
mechanism  ABCD. 
The  verification  for  this  problem  was  accomplished  by  analyzing 
the  mechanism  with  IMP  (Integrated  Mechanisms  Program)  [25],  which  has 
been  implemented  at  the  University  of  Florida  by  the  author.   IMP  is 
primarily  the  product  of  Dr.  John  Uicker,  University  of  Wisconsin-Madison, 
and  performs  kinematic,  static,  dynamic,  time  response,  equilibrium  and 
vibration  analyses  for  virtually  any  type  of  planar  or  spatial  mechan- 
ism via  a  special  user-oriented  language. 

The  dynamic  solutions  are  obtained  by  numerically  integrating 
Hamilton's  canonical  equations  which  are  automatically  generated  from 
the  user's  input  statements  describing  the  mechanical  system. 

IMP  provides  up  to  second  order  analyses  (acceleration)  which 
is  true  of  most  dynamic  analysis  programs.   The  mechanism  has  also  been 
verified  by  the  author's  closed  form,  exact  analysis  program  within 
MECSYN. 

This  analysis  provides  solutions  up  to  the  fourth  order  (jounce) 
to  correspond  with  the  synthesis  capabilities  but  for  reasons  to  be 
disclosed  later  in  this  chapter,  only  performs  analysis  for  the  speci- 
fied precision  conditions.   IMP  was  used  to  provide  a  continuous  anal- 
ysis of  the  mechanism  over  the  range  of  interest.   Figure  (3-2)  illustrates 
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this  analysis.   The  three  curves  are  the  position,  velocity  and  accel- 

(k)  00 

eration  of  the  output  (  $ . (t))  versus  the  position  of  the  input  (  8  (t)). 

Note  that  these  curves  are  not  plotted  versus  time  but  represent  the 
output  of  the  mechanism  when  the  input  has  the  velocity  and  acceler- 
ation prescribed  at  the  first  position.   The  specified  precision  condi- 
tions are  circled  in  Fig.  (3-2). 

Again,  this  illustrates  the  difference  between  geometric 
kinematic  synthesis  and  dynamic  specs  synthesis.   These  curves  would 
be  sufficient  to  describe  the  behavior  of  a  kinematic  function  generator. 

Rather  than  time  derivatives,  they  would  be  if  vs.  9,  di)>/d8  vs.  8  and 

2    2  n 

d  <)>/d8  vs.  8.   To  completely  describe  the  dynamic  specs   function 

generator"  two  sets  of  graphs  must  be  shown,  i.e.,  output  properties  vs. 
time  and  input  properties  vs.  time.   Since  many  different  time  depen- 
dent input  functions  may  be  approximated  by  these  chosen  input  specifi- 
cations, the  analysis  cannot  be  completely  defined  for  this  example. 
Of  course,  in  an  applied  problem,  the  input  specifications  are 
chosen  to  approximate  the  available  time  dependent  input  function,  as 
will  be  demonstrated  in  Chapter  IV.   In  many  cases  the  available  input 
will  be  close  to  constant  speed.   However,  if  the  system  power  shaft 
has  a  known  variation,  dynamic  specs  synthesis  has  the  capability  of 
designing  the  mechanism  for  this  input  as  illustrated  by  these  examples. 

As  was  stated  previously,  this  method  of  four-bar  "function- 
generator"  synthesis  leads  to  a  second  independent  set  of  solutions. 
Step  1  and  Step  2  are  repeated  exactly  as  before.   Instead  of  Step  3, 
proceed  as  follows: 
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(rt  6=45° 


S=l0   9=  is 
<j>=2       <ji=4    $=8 


Figure  3-2  Four-Bar  Output  Link  Properties 
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Step  3' ■   Step  2  resulted  in  a  four-bar,  ABC"D,  whose  coupler  link 

00    00    00 

BC'C"  has  the  dynamic  specs  {  x.(t),   y.(t),   <t.(t)}  , 

00 

while  input,  AB  has  the  specs  8 . (t)  .   One  might  notice 

immediately  from  the  analysis  of  the  path  cognates  for  the 

motion  of  coupler  point  C  in  the  MECSYN  output  that  one 

(k) 
of  the  path  cognates  has  the  output  specs   $.(t),  while  the 

00  J 

input  has  the  specs   6.(t).   This  is  due  to  the  angular 
properties  associated  with  path  cognates  as  described  in 
Chapter  I.   This  path  cognate  mechanism  is  another  solu- 
tion to  the  dynamic  specs  four-bar  "function  generation" 
problem.   The  solution  corresponding  to  the  previously 
synthesized  four-bar,  ABC"D,  is  illustrated  in  Fig.  (3-3) 
along  with  the  "function  generator"  synthesized  in  Step  3 
(ABCD)  and  is  given  by  DFGH. 

D  =  (-1.535,  -1.464)     F  =  (1.659,  0.837) 
G  =  (1.498,  1.032)       H  =  (0.299,  0.016) 
These  two  four-bars  have  only  point  D  in  common  and  since  they 
are  "function  generators"  they  may  be  moved  anywhere  in  the  plane  and 
still  serve  their  purpose.   Therefore,  they  are  independent  solutions. 
For  this  example,  in  Step  2  there  were  four  solutions  for  the  con- 
straint links,  one  of  which  was  the  chosen  link  AB.   Thus,  there  were 
three  possible  solutions  by  angular  cognates  and  three  possible  solu- 
tions by  path  cognates.     As  is  the  case  in  any  mechanism  synthesis 
problem,  especially  applied  synthesis,  many  solutions  must  be  rejected 
for  various  reasons,  among  them,  dimensions,  force  transmission  proper- 
ties, closures,  interference,  physical  compatibility  of  mechanism 
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D=  (-1.535,  -1.464) 
F=(l.659, 0.337) 
G=(l. 498, 1.032) 
H=  (0.299,0.01  6) 


Figure  3-3  Four-Bar  Function  Generator  by  Path  Cognates 
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sections  or  loops,  link  length  ratios,  etc.   For  multi-link  mechanisms, 
these  problems  are  compounded  but  fortunately  more  solutions  are  usually 
available.   It  is  desirable,,  then,  to  have  as  many  solutions  as  possible 
without  giving  up  any  specification  ability.   This  is  especially  true  of 
dynamic  specs  synthesis,  since  it  is  necessary  to  specify  at  least  one 
acceleration  to  treat  dynamics.   Although  many  of  the  six-link  synthesis 
methods  developed  by  the  author  [14]  allow  specification  of  20  conditions, 
they  are  divided  among  two  or  more  links.   Also,  most  of  these  synthesis 
methods  allow  the  choice  of  a  parameter  such  as  a  pivot  location  or  a 
link  length  which  may  be  varied  to  obtain  different  solutions  for  opti- 
mization without  affecting  the  maximum  number  of  dynamic  condi- 
tions specified.   When  this  can  be  done,  as  in  this  four-bar  example, 
it  may  be  desirable  to  reduce  the  specification  set  to  a  four-position 
problem  obtaining  an  infinity  of  solutions.   All  the  synthesis  methods 
are  still  valid.   A  four-position  dynamic  specs  synthesis  curvature 
transform  has  been  included  in  the  latest  version  of  MECSYN  which  uses 
the  same  D   coefficients  and  contains  ordinary  kinematic  synthesis  as 
a  special  case. 

B.   Synthesis  of  a  Stephenson  2  Six-Link  Mechanism 
to  Coordinate  the  Angular  Dynamic  Specs 
of  Three  Links 

A  Stephenson  2  six-link  mechanism  (Fig.  (3-4a))  will  be  synthe- 
sized to  coordinate  the  angular  dynamic  specs  of  the  output  and  a  cou- 
pler with  the  input  of  Case  1  (PPPPP) .   This  means  specification  of  the 
angular  velocity,  acceleration,  jerk  and  jounce  for  three  links  in  one 
location.   Either  9  or  6  may  be  the  input.   In  this  example,  9,  will 
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a. 


A 


Figure  3-4  Inversion  of  Stephenson  2  to  Stephenson  3 
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be  used  as  the  input,  since,  in  the  author's  experience,  it  is  much  more 
difficult  to  find  a  ternary  input  link  (link  5)  which  will  rotate  full 
cycle,  monotonically.   The  length  of  link  3  can  be  arbitrarily  chosen 
thereby  increasing  the  likelihood  that  it  will  satisfy  this  condition. 
That  is  not  the  intent  in  this  example,  but  in  many  instances  of  prac- 
tical significance,  it  may  be  a  requirement.   As  for  the  coupler, 
either  link  2,  1  or  6  may  be  chosen  merely  by  using  the  path  cognate 
feature  of  MECSYN  in  one  design  step.   Since  this  problem  involved  spec- 
ification of  up  to  the  fourth  derivative  with  respect  to  time,  it  will 
be  verified  by  the  author's  closed  form  analysis  within  MECSYN.   This 
analysis  will  be  illustrated  along  with  the  MECSYN  synthesis  steps. 
Step  1.   This  problem  will  be  solved  by  converting  the  specifica- 
tions for  the  mechanism  of  Fig.  (3-4a)  into  equivalent 
specifications  for  one  of  its  inversions,  the  Stephenson  3 

in  Fig.  (3-4b). 

(k)      (k) 
Prescribe  the  angular  dynamic  specs,   {  6,.(t),   6   (5), 
(k)  l2  5J 

63..  (t)},  and  the  location  of  points  A,  B  and  C  as  in 


Fig. 

(3- 

■4a). 

*  J 

k 

(k) 
62j 

(t) 

(k) 
65j 

(t) 

(k) 

(t) 

0  0 

0 

0° 

0" 

90° 

1  0 

1 

6 

2 

5 

rad/sec 

2  0 

2 

7 

4 

10 

rad/sec 

3  0 

3 

2 

8 

15 

rad/sec 

4  0 

4 

3 

8 

18 

rad/sec 

A 

-  (0 

,  0) 

B  =  (1, 

0), 

C  = 

(1,  1). 
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The  actual  initial  angle  of  links  2  and  5  cannot  be  obtained  by 
synthesis  so  they  are  nominally  specified  as  zero. 

Using  the  inversion  feature  of  MECSYN  with  these  specifications, 
the  equivalent  specifications  are  computed  for  point  C  and  link  2  in 
Fig.  (3-4b)  from  the  following  equations: 

00   dk*  . 


*  .  -  9   -  9   +  8,n        where   j,  .  =  51  ,  etc. 

mj    mj    5 j    50  vmj     ,  k 

at 

(k)    (k)    (k)  m  =  2,3,4 

*  .  -   9  .  -   fl  j  =  0 

mJ     mJ     5j  k  -  1,2,3,4. 


(3-2) 


(3-3) 


These  represent  the  inversion  from  the  Stephenson  2  to  the 

Stephenson  3: 

(k)       k 

X.(t)  =  — -  {AB  cos  iji,  +  BC  cos  iji  +  X  } 

J  dtK  J  A  t=t. 

J 

(k)  k 

Y.(t)    =  — r  {AB   sin  i|>,   +  BC   sin  *0  +  Y„ } 
J  dtk  4  3  A|t=t 

i 

which  are  the  equivalent  specifications  for  point  C  in  the  Stephenson  3 
inversion  of  Fig.  (3-4b). 

The  designer  need  not  be  concerned  with  these  equations  since 
MECSYN  carries  out  these  steps  automatically.   The  original  and  the 
equivalent  specifications  are  illustrated  in  Table  3-1  from  MECSYN. 

Step  2.   Using  the  equivalent  inverted  dynamic  specifications, 

(k)      (k)      (k) 
{  X.(t),   Y.(t),    i(in.  (t)},  MECSYN  continues  to  find  two 

constraint  links  (1  and  6)  as  in  Fig.  (3-4b)  with  the 
dynamic  specs  curvature  transform  of  Chapter  II.   The  con- 
verted specifications  and  the  solutions  are  shown  in 
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Table  3-2.   Links  1  and  6  or  EF  and  GH  are  given  by 

E  =  (0.865,  0.810),   F  =  (7.031,  -4.621) 

G  =  (0.698,  -0.683),  H  -  (2.309,  1.191). 

Table  3-3  shows  the  solution  again  for  the  four-bar  EFGHC 

(Fig.  (3-4b))  along  with  its  two  path  cognates,  Grashof's 

mobility  criteria,  link  length  ratios  and  the  analysis  of 

all  three  mechanisms. 

(k)        (k) 
As  was  mentioned  previously,  if   8  . (t)  or  81.(t)  was  specified 
(k)  bJ         lj 

instead  of  8   (t),  one  of  the  path  cognates  illustrated  in  Table  3-3 

would  be  chosen  in  place  of  the  original  mechanism  EFGH  (or  A0,B0„  in 

A   D 

Table  3-3). 

Step  3.   The  Stephenson  3  mechanism  EFGHCBA  synthesized  in  Step  2 
is  inverted  by  releasing  link  FHA  and  fixing  link  AB. 
The  result  is  a  Stephenson  2  mechanism  which  satisfies  the 
desired  specifications.   This  mechanism  is  shown  in 
Fig.  (3-5). 
In  order  to  provide  analysis  verification  for  examples  of  this 
type,  the  author  has  considered  several  types  of  analysis  methods  to  be 
included  as  an  integral  part  of  MECSYN.   Although  many  different  anal- 
ysis methods  are  available,  they  were  not  suitable  for  various  reasons. 
The  desired  analysis  routine  should  have  the  following  properties: 
(a)  in  keeping  with  the  closed  form  nature  of  MECSYN  it  should  be  a 
closed  form  solution  (non-iterative),  (b)  it  should  provide  solutions 
for  up  to  fourth  order  derivatives,  (c)  it  should  be  equally  applicable 
to  all  types  of  mechanisms  capable  of  being  synthesized  by  MECSYN,  and 
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(d)  since  the  complete  input  cycle  is  not  always  known,  it  should 
provide  solutions  for  mechanisms  at  the  precision  points  for  the  order 
of  derivatives  specified  automatically  by  case  number. 

The  routine  developed  to  satisfy  these  conditions  does  not 
overshadow  the  rest  of  the  MECSYN  in  size  and  yet  provides  enough  infor- 
mation about  the  specified  conditions  to  verify  the  mechanism  and  pro- 
vides analysis  of  similar  properties  for  all  other  links  of  the  mechan- 
ism.  The  designer  can  observe  how  other  links  of  the  mechanism  are 
behaving  at  the  precision  locations  in  addition  to  verifying  the  solu- 
tion.  If  additional  full  cycle  information  is  desired,  IMP  or  the  dyad 
package  analysis  [26 ]  may  be  used.   This  analysis  routine  serves  only 
the  needs  of  the  synthesis  process. 

Since  most  analysis  programs  approach  each  different  mechanism 
in  a  different  fashion,  this  type  of  analysis  was  prohibitive  in  size. 
The  routine  to  be  described  here  follows  the  modular  synthesis  approach 
of  MECSYN. 

Figure  (3-6)  shows  the  path  cognates  of  an  original  set  of  two 
dyads,  CAO  and  CBO  which  are  partial  results  of  a  synthesis  problem. 
The  remaining  four  dyads  represent  the  first  and  second  cognate  solu- 
tions of  MECSYN  (Table  3-3) .   The  locations  of  all  these  points  for  the 
first  specified  position  are  known  after  a  synthesis  and  cognate  calcu- 
lation step  of  MECSYN.   Depending  on  the  case  number,  a  series  of 
vectors  for  position,  velocity,  acceleration,  jerk  or  jounce  is  applied 
to  point  C  which  consists  of  all  the  path  specifications  for  point  C  in 
the  synthesis  step  (Fig.  (3-6)).   For  the  synthesis  step  to  be  valid, 
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Figure  3-6  Analysis  of  a  Four-Bar  and  Its  Path  Cognates 
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coupler  links  formed  by  AC  and  BC  must  have  Identical  angular  properties 

00 
equal  to  the  angular  properties   #.(t)   specified  for  the  coupler  link 

(k)     (k) 
when  point  C  has  the  specified  state   {  x.(t),   y.(t)}.   If  this  is  the 

case,  the  two  dyads  form  a  valid  four-bar  solution  for  the  synthesis 
problem  and  dimension  AB  is  constant.   Also,  note  each  cognate  dyad  con- 
straint link  D0d  and  FO  must  have  precisely  the  same  angular  prop- 
erties as  those  specified  for  the  coupler.   These  corresponding  proper- 
ties are  enclosed  in  rectangles  in  Table  3-2  and  Table  3-3.   Other 
similarities  between  the  angular  properties  of  links  in  path  cognates 
as  described  in  Chapter  I  may  also  be  observed  in  Table  3-3. 

Referring  to  Fig.  (3-7),  the  analysis  is  formulated  as  follows: 
The  location  of  the  points  A  and  0  in  the  first  position  for  all  the 

dyads  of  Fig.  (3-6)  is  determined  by  the  synthesis  and  path  cognate 

(k)      (k) 
steps.   The  dynamic  specifications       {  x  . (t) ,   y  .(t)}   for  the 

1      •    n  ^     (k)CJ 

coupler  point  C  are  represented  in  Fig.  (3-7)  by   C . (X  ,Y  ).   Thus  the 

(k)      (k)        J   2  2 

known  quantities  are  A,  0,  and  {  x  .(t),   y   (t)}. 
A        cj      ■'cj 

^X1'Y1^  are  the  coordinates  of  0  and  (X  ,Y  )  are  the  coordinates 
of  the  first  position  of  C.   The  desired  properties  are 

00         (k) 
*y(t)  and  *   v  z' 

for  both  closures  of  each  dyad  (whether  point  A  is  above  or  below 
line  0AC) .   Both  closures  are  necessary  since  frequently  one  or  more 
dyads  may  change  into  another  closure  while  satisfying  the  prescribed 
conditions.   This  does  not  mean  that  the  mechanism  will  have  closure 
difficulties  in  satisfying  the  specifications.   For  example,  the  four- 
bar  0AACB0B  of  Fig.  (3-6)  will  have  closure  difficulties  if  dyad  ABO 
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Figure  3-7     Dyad  Analysis 
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must  go  into  a  different  closure  when  AO  is  the  input,  regardless  of 
what  closures  the  dyads  CAO.  and  CB0_  are  going  into. 

A  n 

This  closure  problem  is  more  difficult  to  define  for  six-bars 
in  terms  of  dyad  closures.   It  appears  that  a  test  for  closures  must 
be  written  independently  for  each  mechanism  rather  than  relying  on  a 
modular  approach  as  is  done  for  all  other  components  of  MECSYN. 
Solution  of  this  closure  problem  is  necessary  to  provide  complete  com- 
puter selection  of  optimum  mechanisms.   If  another  analysis  routine  is 
used  such  as  IMP  or  the  dyad  package,  it  will  be  immediately  obvious 
whether  a  mechanism  satisfies  all  prescribed  conditions  in  one  closure 
since  these  analyses  are  performed  for  each  mechanism  separately  and  in 
the  original  closure.   This  approach  would  be  prohibitive,  costwise,  for 
evaluation  of  a  large  number  of  synthesized  mechanisms  to  obtain  optimum 
solutions. 

The  angles  i|>-  and  i|i  in  Fig.  (3-7)  are  given  in  the  finitely 
separated  positions  by: 


where  ±  det 

ermines  the  closure 

and 

0  <  3 

<  n 

o  <  s2 

<  n 

*2j    "  «J    *   hi 


D2   -    (X2.-Xl)2  +    (Y^-y/ 

Bj  -   cos"1    [(L2   -   L2  +  D2)/2L1D]  h1   =    |A0A| 

-1  2  2  2  L?   =    IACI 

B     =   cos        [(i/  -  L,   -  DV2L„D] 

12  D     -    |0AC| 


(3-4) 
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The  vectors  ty_  and  r_  and  the  matrix  A  are  defined  by: 


Vtr 

'x2.(t)  -x1 

*2j(o 

r  = 

Y2.(t)    -Yx 

-L     sin  i|i      (t) 

L2   sin  <i2At) 

1^   cos   ^x    (t> 

-L     cos   i|i_.  (t) 

th 


All  the  zero   order  properties,  iji.(t),  i|i„  (t),  are  given  by 
equation  (3-4).   The  higher  order  properties, 


J\At) 


d  *2,(t) 
and  *s 


dt  dt 

are  determined  from  the  following  equations  (where  j^  represents  dj|j/dt, 
etc.) : 

1  =  a-1  i 


i  =  A"    [t  -  A  ii 


„-i 


2  =  a     [r*  -  2  A  $  -  A 1] 


(3-5) 


-lr  " 


±     =A   [  r-  3Ai-  3A|-A  i] 


A  representative  example  (case  1)  of  this  analysis  is  shown 
(k) 
in  Table  3-3.   The  properties  \ji    .  (t)   correspond  to  the  columns  labeled 

J   00 
constraint  link  and  the  properties   i/>_  .  (t)  to  those  columns  labeled 

coupler  link.   The  columns  of  the  analyses  correspond  to  the  mechan- 
isms listed  directly  above  them.   More  detailed  information  about  the 
use  of  this  analysis  routine  may  be  found  in  the  description  of  MECSYN 
to  follow  (under  the  heading  subroutine  ANALYS). 
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For  the  previous  Stephenson  2  example,  Tables  3-2  and  3-3 
give  verification  of  the  inverted  four-bar  which  satisfies  the  inverted 
specifications  for  point  C.   What  remains  to  be  verified  are  the 
inverted  specifications  for  point  C,  of  the  dyad  ABC  in  Fig.  (3-4b) . 

However,  this  can  be  immediately  verified  as  was  done  in 

Section  A  for  the  four-bar  "function  generator."  That  is,  choosing 

arbitrarily  a  dyad  ABC  and  the  angular  specifications  desired  for  this 

GO       (k) 
dyad,  i|i   (t),   *31(t)  as  in  Fig.  (3-4b) .   MECSYN  is  used  to  find 

the  equivalent  specifications  for  point  C.   Synthesizing  constraint  links 

(k)      (k)      (k) 
for  this  dyad  using  the  specifications  {  X  .(t),   y   (t),   ij;   (t) } 

CJ       cj       3j 

gives  two  or  four  constraint  links  in  general.   For  the  computed  speci- 
fications of  point  C  to  be  correct,  one  of  these  solutions  must  be  the 
arbitrarily  chosen  constraint  link  AB.   In  this  fashion  the  routine  that 
computes  the  inverted  specifications  of  point  C  has  been  verified  for 
all  seven  cases. 

Thus  there  are  three  independent  self-checking  routines  within 
MECSYN;  the  routine  which  computes  the  inverted  specifications  for 
point  C  of  a  chosen  dyad  ABC,  the  dynamic  specs  synthesis  curvature 
transform  and  the  analysis  routine  described  previously  in  this  chapter 
which  computes  angular  properties  given  the  path  specifications.   Since 
these  verification  techniques  are  applied  modularly  just  as  the  syn- 
thesis procedure  is  applied  to  a  mechanism,  any  mechanism  that  can  be 
synthesized  can  be  immediately  verified.   An  additional  benefit  arises 
from  this;  i.e.,  if  any  errors  are  made  in  specification  during  the 
design  steps,  the  errors  are  quickly  revealed  by  comparing  these  self- 
checking  features. 
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C.  Dynamic  Specs  Synthesis  of  a  Geared 
Five-Link  Mechanism 

To  illustrate  the  versatility  of  MECSYN  a  geared  five-link 
mechanism  will  be  designed  for  prescribed  dynamic  specs  of  coupler 
link  2  with  input  link  5  for  Case  4  (PPP-P-P)  as  illustrated  in 
Fig.  (3-8a).   This  implies  specification  of  velocities  and  acceler- 
ations at  the  first  position  and  two  additional  finite  positions. 
Some  preliminary  formulations  are  necessary  for  the  design  of  the  geared 
five-bar.   The  gear  ratios  as  shown  in  Fig.  (3-8a)  are: 

M  =  r  /r.  ~N  M  =  -r  /r.  > 

a  b  I  a  b  I 

>   with  idlers  >   without  idlers 

N  "  VrJ  N  =  -VrJ 

Q  =  MN  =  gear  ratio  between  points  A  and  C. 
The  angular  relationships  between  links  2,  3  and  gear  B  are 
A62  =  N(A0B  -  A93)  +  A93 

A8_  =  M  A8_ 

A8   -  N(M  A9   -  AS  )  +  A6 

A92  =  Q  A65  +  (1-N)  A63 

The  higher  order  angular  properties  are  obtained  by 


^k  (92i-e20)  =  ^k 
dtfc   ZJ  ZV  dt 


e2  =  q  e5  +  (i-N)e3 


92  =  Q  65  +  (1-N)93 


Q(85.  -  e50)+(l-N)(83-930) 
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Figure  3-8  Geared  Five-Bar  Inversic 
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Thus  the  relation  of  the  input  link  5  to  links  2  and  3  is  given  by 
A65   =  Q  U92  "  (1-N)ae3} 
§5  =  Q   {52  "  (1-N)83}  (3_6) 

e5  =  ^  {§2  -  d-N)e3}  . 

A  constant  speed  input  will  be  selected,  i.e.,  B  =  0,  8  =  0.1  rad/sec. 
Two  free  parameters  are  available  which  allow  solution  of  this  problem 
for  dynamic  state  specifications,  namely,  the  gear  ratios  M  and  N. 
In  the  case  of  geometric  kinematic  synthesis,  the  specified  motion  of 
link  2  would  be  achieved,  regardless  of  what  input,  8   is  supplied. 
If,  however,  it  is  desired  that  coupler  link  2  satisfy  a  prescribed 
inertia  torque  and  force,  a  specific  velocity  and  acceleration  must  be 
supplied  at  the  input,  9  .   The  necessary  input  may  be  found  by  analysis 
after  synthesis  and  provided,  but  it  is  not  usually  feasible  to  provide 
an  arbitrary  acceleration  and  velocity  to  the  input.   Instead  the  mechan- 
ism should  be  synthesized  to  be  driven  with  the  available  velocity  and 
acceleration  of  the  power  shaft.   The  parameters  M  and  N  are  determined 
in  the  synthesis  procedure  to  allow  the  desired  input  velocity  and 
acceleration. 

N  is  chosen  to  give  the  input  acceleration  condition  at  position 
zero,  6  =  0. 

0  =  |  {S2  -  (1-N)S3) 

n  =  l  -  (e'2/e'3)  .  (3_7) 
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Q  is  chosen  to  coordinate  the  desired  input  velocity  6.  =  0.1  rad/sec, 
with  the  specifications  of  link  2. 

95=i(92-  (1-N)93} 

Q  =  — (e,  -  (i-N)e,} 

85  L  i 

q  =  io  {S2  -  (e'2/e'3)93}  .  (3-8) 

(k)  (k) 

The  properties   9   (t)  are  prescribed  and  the  properties  9   (t) 

arise  from  the  analysis  routine  in  MECSYN  after  the  first  synthesis 

step.   N  and  Q  can  then  be  found. 

Step  1.   Prescribe  the  dynamic  specs  of  link  2, 

(k)       (k)       (k) 
{  Xp.Ct),   YD-(t).    92.(t)},  tht  velocity  9   and  the 

acceleration  8  of  the  input  and  the  location  of  point  A 


(Fig- 

(3- 

•8a)). 
(k) 

(k) 

(k) 

*■    i 

k 

v° 

YDj(t) 

e2(t) 

0  0 

0 

2  in. 

2  in. 

0° 

1  0 

1 

2  in. /sec 

2  in. /sec 

4  rad/sec 

2  0 

2 

1  in. /sec 

-3  in./sec 

6  rad/sec 

3  1 

0 

2.5  in 

2.3  in. 

10° 

4  2 

0 

3  in. 

2.5  in. 

15° 

95=0 

1  rad/ 

sec 

2       A  =  (0,  3). 
8=0  rad/sec 
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Using  the  dynamic  specs  curvature  transform  of  MECSYN,  a  suitable  con- 
straint link  BC  is  found  and  shown  in  Fig.  (3-9) 


B  =  (1.248,  3.076) 
(k) 


C  =  (0.755,  4.200). 


The  values  of  8  .  (t)  for  this  link  are  obtained  from  MECSYN's  analysis. 

Step  2.   The  synthesis  is  carried  out  similar  to  that  of  a 

Stephenson  2  mechanism.   Referring  to  Fig.  (3-8),  link  5 

will  be  fixed  and  link  4  released.   The  inverted  properties 

(k)       (k) 
of  dyad  ABC   (  <j> ,    (t) ,   4>   Ct> )  will  be  found  from 


(k) 

e5j(t), 


(k) 


(k) 


3,.(t)   and   6,.(t).   Using  the  inverted 


00 


00 


(k) 


dynamic  specs  {  ^,,(t),   Yc.(t),   ^.(t)}, 

a  constraint  link  EF  will  be  found  as  shown  in  Fig.  (3-8b) . 
For  the  geared  five-bar  some  intermediate  calculations  must  be  carried 
out  first. 

From  the  MECSYN  analysis 


S.      j  k 


GO 

e2.(t) 


(k) 


0  0  0  119.49271° 
10  1  4.0  rad/sec 

2  0  2  6.0  rad/sec 

3  10  129.49271° 

4  2  0  134.49271° 
Using  Eqs.(3-7)  and  (3-8)  and  the  above,  N  and  Q  can  be  calculated  as 

N  =  0.34598,  Q  =  0.43408,    M  =  Q/N  =  1.22574. 


1 

>3j(t) 

113 

.64762 

6 

.05118 

rad/sec 

9 

.17404 

rad/sec 

106, 

.83605° 

90, 

.90280° 

Specifications 


76 


pos. 

der. 

xD 

Yd 

92 

e5 

0 

0 

2 

2 

0° 

0 

1 

2 

2 

4 

0.1 

0 

2 

1 

-3 

6 

0 

1 

0 

2.5 

2.3 

10" 

2 

0 

3.0 

2.5 

15° 

A=(0,3) 


I- 


o 

I 


o 

2 


M=r0/rb=  1.2  2  S 
N  =rb/rc=  0.34  6 


A=  (0,3) 

B  =(1.248,3.073) 

C  =(0.755,4.200) 

D=(2,2) 

E  =  (I.8I6,2.523) 

F  =(-0.307,2.2  49) 


Figure  3-9  Synthesis  of  Geared  Five-Bar  -  Position  0 
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Figure  3-10  Synthesis  of  Geared  Five-Bar  -  Positions  1,  2 
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Applying  Eq.  (3-6)  will  give  the  angles,  8 

.    i65  "  85j  "  950  "  i   [A92  "  Cl-H)463] 


(951  "  950}  =  i951  =  34.08532° 


(952  "  950)  =  A952=  70-447950 


Using  these  specifications 

00 


i 

i 

k 

e5.(t) 

0 

0 

0 

0° 

1 

0 

1 

0.1  rad/sec 

2 

0 

2 

2 
0     rad/sec 

3 

1 

0 

34.08532° 

4 

2 

0 

70.44795° 

(k)      (k) 

along  with  the  previously  listed  specifications  for  8,.  and  8„ 

2j       3j 

with  MECSYN,  the  inverted  specification  set  {  X   (t)  ,  y^..(t),     i|i, .  (t)  } 

Cj       Cj  2j 

is  found. 

(k) 
The  angular  properties  \ji    .  (t)  of  Fig.  (3-8b)  are  related  to 
(k)  raJ 


8 (t)  of  Fig.  (3-8a)  by 


*  .= 

mj 

8  .  -  8..  +  ec„ 

mj           5j           50 

(k) 

(k)           (k) 

=      8    .    -      6 
mj             5j 

=  1... 


for  k  >  0. 


S_tep_3.   Using  the  inverted  specification  set  for  the  dynamic  specs 
of  link  2  in  Step  2,  a  constraint  link  EF  is  synthesized 
with  MECSYN. 


79 


I 

J 

k 

(k) 

xc.Ct) 

(k) 

(k) 
Cj(t) 

0 

0 

0 

0.7554 

in. 

4.2005° 

119.4927° 

1 

0 

1 

-6.6820 

in. /sec 

-3.0540  in 

/sec 

3.9  rad/sec 

2 

0 

2 

7.1076 

in. /sec 

-44.3275  in 

,   2 
/sec 

2 
6.0  rad/sec 

3 

1 

0 

1.4400 

in. 

3.5360  in 

95.4070° 

4 

2 

0 

1.6393 

in. 

2.2788  in 

64.0448° 

E  =  (1.816,  2.523)      F  =  (-0.307,  2.249). 

Step  4.   Reinvert  Che  mechanism  by  fixing  link  4  and  releasing  link  5. 
Applying  the  gearing  M,  N  from  point  A  to  C  gives  the 
desired  geared  five-bar  mechanism  ABCDEF  shown  in  Figs.  (3-9) 
and  (3-10). 

The  following  chapter  presents  solutions  to  two  practical  prob- 
lems using  the  methods  of  Chapter  II  and  this  chapter. 


CHAPTER  IV 


EXTENDED  PRACTICAL  APPLICATION 
OF  DYNAMIC  SPECS  SYNTHESIS 


The  theory  of  Chapter  II  is  applied  in  this  chapter  to  the 
solution  of  two  problems.   First,  the  dynamic  specs  of  a  link  moving 
with  general  coplanar  motion  and  the  corresponding  dynamic  input  specs 
are  prescribed  to  synthesize  a  frequently  needed  transfer  device 
between  linear  and  circular  moving  conveyors.   Both  a  Stephenson  3 
and  a  Stephenson  2  six-link  mechanism  are  synthesized  for  the  same 
transfer  specifications.   Second,  a  Stephenson  3  six-link  mechanism  is 
synthesized  for  torch  or  bandsaw  cutting  of  large,  low  pitch  (diame- 
tral pitch  -  3/8)  gear  teeth.   The  geometric  specs  of  a  coupler  link 
(path  and  angular  properties)  are  coordinated  with  the  corresponding 
geometric  input . 

Today  the  designer  is  faced  with  a  much  wider  range  of  systems 
capable  of  solving  any  given  mechanical  motion  problem.   One  end  of 
this  spectrum  is  represented  by  automatic  control  devices  composed  of 
electronic,  pneumatic,  or  hydraulic  components.   The  other  end  is 
dominated  by  single  degree  of  freedom  systems  such  as  cams  and  linkages. 
The  spectrum  has  been  dramatically  enlarged  during  the  past  three 
decades  primarily  due  to  the  increased  activity  in  the  design  and  appli- 
cation of  automatic  control  devices.   This  led  to  a  relative  de-emphasis 
of  study  and  research  in  the  mechanism  field  limiting  the  flexibility 
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with  which  the  designer  can  meet  any  given  problem.   The  objective  of 
most  synthesis  procedures  is  to  re-establish  mechanisms  as  an  effective 
and  viable  alternative  to  systems  based  solely  on  control  type  components 

Mechanisms  are  generally  considered  devices  with  limited  adjust- 
ability capable  of  performing  a  given  collection  of  functions  with  high 
precision  and  reliability.   Unless  internal  stresses  and  energy  stor- 
age lead  to  vibration  and  consequent  destruction,  these  systems  can  be 
operated  at  high  speeds,  can  drive  large  masses,  and  overcome  large 
forces.   The  cam  is  the  most  flexible  information  storage  device  among 
mechanisms  but  it  is  the  least  capable  of  fulfilling  the  above  require- 
ments because  of  its  structural  weakness.   Linkages  can  be  made  to 
handle  very  high  loads  but  they  can  only  meet  a  limited  range  of  motion 
functions.   Consequently,  designers  of  linkages  must  understand  the 
requirements  of  their  output  functions  to  a  greater  degree  and  be  very 
careful  not  to  overspecify  the  problem  because  the  linkage  has  a  limited 
number  of  free  design  parameters. 

These  limitations  of  linkages  have  prevented  their  broad  use 
in  recent  years.   Design  procedures  as  presented  here  are  intended  to 
strengthen  the  designer's  capability  in  creating  a  mechanism  to  satisfy 
a  needed  function  or  group  of  functions.   Utilizing  dynamic  specs  syn- 
thesis linkages  can  now  be  asked  to  perform  combinations  of  functional 
requirements  such  as  position,  velocity,  and  acceleration  of  a  mass  or 
position,  energy  and  power  associated  with  a  work  function.   Further, 
recent  developments  in  dynamic  balancing  [27]  make  it  fully  possible 
to  place  complex  linkages  in  a  high  speed  system  and  still  ensure 
smooth  operation.   Springs,  dashpots,  and  masses  can  all  be  added  to 
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the  system  to  enhance  this  smoothness.  Linkages  will  continue  to  con- 
sume more  space  than  other  systems  and  not  all  their  properties  can  be 
controlled  by  direct  synthesis  techniques.  Hence,  the  role  of  optimi- 
zation will  always  prove  dominant  for  critical  design  applications  and 
therefore,  a  full  computational  capability  based  on  synthesis  to  reduce 
the  number  of  free  design  parameters  will  be  required. 

A.   Dynamic  Specs  Synthesis  of  a  Stephenson  3  Mechanism 

as  a  Transfer  Device  Between  Linear  and  Rotary  Conveyors 

A  common  problem  in  materials  handling  is  that  of   transfer- 
ring an  object  from  a  continuously  moving  linear  conveyor  belt  to  a 
continuously  rotating  turret.   If  the  linear  surface  speeds  of  the  two 
devices  are  close  in  magnitude,  the  objects  frequently  can  be  guided 
by  fixed  plows  onto  the  turret  surface.   It  may  be  possible  to  design 
a  mechanism  by  geometric  motion  specifications  to  match  the  required 
path  properties  (i.e.,  tangent,  curvature,  etc.)  should  the  plow  con- 
cept prove  inadequate.   If,  however,  the  linear  surface  speeds  differ 
markedly,  then  it  is  unlikely  that  the  required  matching  of  velocity, 
acceleration,  jerk,  etc.,  would  result  and  dynamic  shock  will  be  expe- 
rienced by  the  object  during  the  transfer  operation.   In  some  cases 
where  the  mass  and  rotational  inertia  of  the  object  is  low  or  if  the 
object  or  machine  function  is  largely  insensitive  to  shock,  this  mis- 
match of  dynamic  specs  may  not  be  deleterious.   On  the  other  hand, 
matching  of  the  dynamic  specs   could  improve  system  precision  of  oper- 
ation or  allow  a  significant  increase  in  production  rates. 

In  this  section,  a  mechanism  will  be  designed  using  dynamic 
specs  synthesis  which  will  match  the  dynamic  specs  of  the  object  to  the 
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dynamic  conditions  it  would  experience  on  the  turret  surface.   Also,  the 
angular  and  linear  velocity  of  the  object  at  pick-up  will  match  the 
velocity  condition  of  the  belt.   Finally,  these  output  dynamic  specs 
will  be  coordinated  with  the  available  input  from  the  power  shaft  of  the 
system. 

The  objects  are  spaced  1  foot  apart  on  the  conveyor  belt  moving 
with  constant  velocity  of  1  foot/sec  with  an  orientation  of  0°  with  the 
horizontal  and  zero  angular  velocity.   The  turret  rotates  at  a  constant 
speed  of  0.5  revolutions/sec  with  an  effective  operating  radius  of  4/tt 
feet.   The  location  of  the  turret  gripper  fingers  at  transfer  is  135° 
from  the  X  axis.  At  this  location  on  the  turret,  the  object  must 
undergo  a  velocity  of  4  ft/sec  with  an  orientation  of  45°  with  the 

X  axis,  an  angular  velocity  matching  that  of  the  turret  or  0.5  revolu- 

2 
tions/sec  (irrad/sec),  and  an  acceleration  of  4tt  ft/sec  directed 

toward  the  center  of  the  turret.   The  input  shaft  has  a  constant  angu- 
lar velocity  of  2tt  rad/sec.   An  object  must  be  transferred  during  each 
1-sec  interval,  while  the  actual  transfer  process  is  to  consume  1/2  of 
this  interval  or  0.5  sec  leaving  0.5  sec  for  the  mechanism  to  return  to 
the  beginning  of  the  cycle.   Two  stages  are  provided  on  the  turret  180° 
apart  to  allow  two  objects  to  be  placed  on  the  turret.   After  the  neces- 
sary turret  operation  on  the  object,  it  may  be  removed  by  an  identical 
mechanism  on  the  opposite  side.   The  component  dynamic  specifications 
for  this  problem  are  illustrated  in  Fig.  (4-1). 

The  stated  synthesis  problem  is  now  presented  in  terms  of  the 
three  design  steps  outlined  in  Chapter  II. 
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{ [3>      t^VA=' 


conveyor 


V  l_ 

0 
A  |_ 


A-  PICKUP 
X  =  1,8  ft 
Y  =-1.0  ft 

$=  0  rod 
8  =Trod 
X  =   l.0ft/4«c 
Y=0  tt/stc 
<f>  =  0  rad/itc 
g  =  27T  red/sec 


B -  RELEASE 

X  =  2.6  « 

Y=  0.9ft 

<#>=7r/4rad 

8=0  rod 

X  =  2.8284  rt/s.c 

Y=2.82B4  ft/j«c 

i£  =  it  rod/^BC 

e=2^rod/3<c 
X  =    8.8858  ft/s.c1 
Y  =-8.8  858  ft/j.c1 
(£  =  0  rad/s»ca 
fl=0  rad/j.c-1 


Figure  4-1     Problem  Specifications 
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Figure  4-2  Transfer  Mechanism  -  Drop-Off  Position 
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Figure  4-3  Transfer  Mechanism  -  Pick-Up  Position 
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Step  1.   Prescribe  the  dynamic  specs  as  given  in  Fig.  (4-1): 


I 

J 

k 

00 

00 

00 

0 

0 

0 

2.6  ft 

0.9  ft 

45  deg 

1 

0 

1 

2.8284  ft/sec 

2.8284  ft/sec 

¥ 

rad/sec 

2 

0 

2 

8.8858  ft/sec2 

-8.8858  ft/sec2 

0 

rad/sec 

3 

1 

0 

1.5  ft 

-  1  ft 

0  deg 

4 

1 

1 

1   ft/sec 

0  ft/sec 

0 

rad/sec 

Using  the  formulation  of  Chapter  II,  consider  Case  3  (PPP-PP) 

and  determine  by  the  dynamic  curvature  transform  a  single  constraint 

crank  BO   to  satisfy  these  specifications.   Fig.  (4-2)  shows  the  final 

choice  for  BO  in  the  first  position  of  the  system.   Note  that  the 

first  finite  position  (j=0)  is  chosen  at  the  drop-off  location  to  agree 

with  the  computer  program  MECSYN  formulation  which  requires  that  any 

finite  position  for  which  k  >  2  be  specified  first.   This  can  always 

be  done  without  loss  of  generality  (even  for  such  cases  as  PPP-P-P) 

since  the  sequence  in  which  the  mechanism  passes  through  finitely 

separated  positions  (FSP)  is  independent  of  the  labeling  required  for 

the  analytical  formulation.   The  numerical  coordinates  of  the  selected 

dyad  CB0  are: 

C  (2.6,  0.9),     B  (0.361,  2.252),     0.  (0.873,  2.079). 

b 
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fA9=  180°  1 j 

\At  =  l/2««c.J 


Figure  4-4  Transfer  Mechanism  -  Path  Properties 


89 


<t>       140 

CdegreeO 


(rac^ec) 
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Figure  4-5     Transfer  Mechanism  -  Angular  Properties 
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Step  2.   Prescribe  the  specs  of  point  C  and  the  input  specs  of  the 
power  shaft  for  the  dynamic  specs  of  a  moving  plane: 


00 

(k) 

00 

I 

J 

k 

xj 

7i 

9. 
J 

0 

0 

0 

2.6    ft 

0.9   ft 

0  deg 

1 

0 

1 

2.8284  ft/sec 

2.8284  ft/sec 

2ir  rad/sec 

2 

0 

2 

8.8858  ft/sec2 

-8.8858  ft/sec2 

2 
0  rad/sec 

3 

1 

0 

1.5   ft 

-1.0  ft 

180   deg 

4 

1 

1 

1.0  ft/sec 

0        ft/sec 

2tt  rad/sec 

Using  MECSYN  find  two  constraint  links  forming  a  four-bar  linkage  that 
satisfy  these  dynamic  specifications.  Path  cognates  of  this  four- 
bar  linkage  will  have  one  crank  (F0f)  which  undergoes  the  desired  input 

(k)  t  (k)   (k) 

specs  9.  while  point  C  experiences  the  specified  specs  {x.,   y .} .  The 

final  choice  of  the  four-bar  cognate  to  satisfy  this  step  is  shown  in 

Fig.  (4-2)  as  GO  ,  FO  ,  C  with  the  numerical  coordinates: 

G  (1.591,  0.395)      0   (0.356,  1.619) 
g 

F  (1.299,  -0.414),    0f  (0.846,  0.033). 

Step  3.   The  resulting  link  systems  of  Steps  1  and  2  have  in  common 
00     (k) 
the  specs  (x.(t),   y.(t)}at  point  C.   Joining  the  dyad  to  the 

cognate  four-bar  gives  the  desired  Stephenson  3.   The  motion 

of  the  mechanism  in  its  two  distinct  design  positions  j=0,l 

is  illustrated  in  Figs.  (4-2)  and  (4-3)  with  accompanying 

point  path  for  C.   Fig.  (4-4)  shows  the  x,y  components  of 

velocity  and  acceleration  for  point  C  over  the  whole  cycle. 

Figure  (4-5)  illustrates  the  angular  position,  angular 
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00 

velocity,  and  angular  acceleration  of  link  CB  as   4  .   All  are 

J 
plotted  versus  the  input  angle  a.   All  prescribed  data  are  indi- 
cated by  circles  on  these  graphs.   The  conditions  as  specified 
were  confirmed  by  two  independent  analysis  procedures  to  further 
validate  the  results.   Note  the  graphs  do  indicate  relatively 
high  local  accelerations  outside  the  control  of  the  synthesis 
process.  As  mentioned  in  the  introduction,  such  phenomena  can 
be  minimized  only  by  careful  optimization. 

B-   Dynamic  Specs  Synthesis  of  a  Stephenson  2 
Mechanism  as  a  Transfer  Device  Between 
Linear  and  Rotary  Conveyors 

Using  the  same  problem  description  and  data  given  for  the 
Stephenson  3  problem,  a  Stephenson  2  mechanism  will  be  synthesized  to 
illustrate  the  application  of  the  four  tools  of  synthesis  in  a  modular 
fashion  to  a  more  difficult  mechanism  synthesis.   The  first  objective 
in  the  synthesis  procedure  is  to  convert  the  desired  motion  specifica- 
tions into  numerical  specifications  suitable  for  use  in  the  dynamic 
curvature  transform  central  to  MECSYN.   The  specifications  for  the 
Stephenson  3 

/oo  00     00     00   "\ 

lxcj(t)'   V0,   Vc)'   9j(t)J 

are  now  described  by  a  set  of  symbols  corresponding  to  those  shown  in 
Fig.  (4-6a)  where  point  D  replaces  C  as  the  path  generator  with  the 
specifications 

foo  00     00     00  1 

IV0'    Vt}-    e5j(t)'    V°/  • 
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V4 


U 


Figure  4-6  Two  Inversions  of  the  Stephenson  Six-Link  Chain 
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In  addition  the  location  (U.,V  )  of  pivot  A  can  be  prescribed  in  the 
reference  £(U,V). 


Step  1.   The  first  requirement  is  to  determine  the  link  configuration 
BCD  as  a  segment  of  links  2  and  3  of  Stephenson  2.   This  link 
configuration  is  identical  to  0  BC  in  the  Stephenson  3  so  that 
the  means  to  satisfy  this  step  is  the  same  as  Step  1  for  the 
Stephenson  3.   It  is  not  necessary  that  the  same  link  solution 
from  the  dynamic  curvature  transform  be  used  as  before.   The 
same  solution  is  used  here,  however,  to  illustrate  the  similar- 
ity between  the  two  mechanism  problems.   The  numerical  coordi- 
nates of  this  solution  are  given  as: 

A  (1.25,  1.25),     D  (2.6,  0.9) 

C  (0.361,  2.252),    B  (0.873,  2.079)  . 

Step  2.   In  this  part  of  the  synthesis,  the  Stephenson  2  is  inverted  by 
releasing  link  4  and  fixing  link  5  giving  a  Stephenson  3  as 
illustrated  in  Fig.  (4-6b).   It  follows  from  the  analytical 
rules  of  inversion  that  the  specs  of  point  C  and  the  angular 
specs  of  link  2  can  be  obtained.   This  means  that  the  dynamic 
specs  of  link  2  of  the  Stephenson  3  can  be  prescribed.  The  chosen 
point  A  along  with  the  synthesized  points  B  and  C  from  Step  1 
form  the  dyad  structure  ABC.   This  dyad  can  be  analyzed  to  pro- 
vide the  specs  of  point  C  after  the  inverted  angular  specs 
(k)    (k) 
iji,  .,  ii..    are  determined.   MECSYN  provides  this  dyad  analysis 

and  angular  inversion  for  dynamic  synthesis  according  to  case 

number  of  DMSP.   The  angular  inversion  {6  ■*  i|>  }  between  the 

mj   mj 
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Stephenson  2  and  the  Stephenson  3  for  link  m  for  position  j  is 
given  by  the  formula 

♦  .  =  e  .  -  e_.  +  ec„ 

mj    mj     5j    50 

where  e50  is  the  initial  value  of  65  in  the  invariant  reference  system 
E(U,V).   The  derivative  of  this  displacement  relation  gives  the  angular 

00   (k) 

inversion  of  specs  {  9  .  •+■  ill   }  as 
mj    rmj 

(k)     (k)     (k)     (k) 

**j  =  9mj  "  S5j  +  V    m  -  2>3>4  • 

Note  that  all  angles  are  measured  from  the  same   (U,V)  reference  system 

for  both  inversions.   The  inverse  of  the  above  process  {  ill   +  6   }  is 

mj     mj 

(k)     (k)     (k)     (k) 

8  .  =  <L    .    -   ill,  .  +  ill, 
mj    Tmj    r4j    ^40 

since  link  4  is  fixed  in  the  Stephenson  2.   Now  we  may  obtain  by  direct 

analysis  the  specs  of  link  C  for  the  Stephenson  3  as 

(k)     dk   _         _ 
XDJ  =  — £  {AB  cos  i|i4  +  BC  cos  iji  +  UA>  , 


dt 


1  t=t. 
J 


(k)     dk  _         _ 

yDj  =  — k  {AB  sin  *4  +  Bc  sln  *3  +  Vi 

dt  lt=ti 

J 

This   completes   the  desired    dynamic    specs 

00    (k)    (k) 

*Dj  '  yDj  '  *2i '  f or  the  lnverted  properties  of  link  2  in 
Fig.  (4-6b). 
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Step  3.   Using  the  same  numerical  specifications  from  Section  1,  the 
initial  input  data  are 

(k)  (k) 


£ 

J 

k 

e5j(t) 

e2.(t) 

0 

0 

0 

0  deg 

45  deg 

1 

0 

1 

2ir  rad/sec 

tt  rad/sec 

2 

0 

2 

2 
0  rad/sec 

2 
0  rad/sec 

3 

1 

0 

270  deg 

0  deg 

4 

1 

1 

2tt  rad/sec 

0  rad/sec 

with  the  exception  that  0   was  changed  from  180c  to  270°  to  achieve 
a  somewhat  better  solution.   By  the  formulation  of  this  step,  the  dynamic 
specifications  for  link  2  in  the  Stephenson  3  are: 


00 

(k) 

(k) 

l 

J 

k 

v° 

v° 

*2J(t) 

0 

0 

0 

0.3612  ft 

2.2521  ft 

45  deg 

1 

0 

1 

4.8770  ft/sec 

1.3796  ft/sec 

-it  rad/sec 

2 

0 

2 

13.2291  ft/sec2 

-43.9557  ft/sec2 

2 

0  rad/sec 

3 

1 

0 

0.9608  ft 

0.8730  ft 

-270  deg 

4 

1 

1 

-2.3689  ft/sec 

2.8169  ft/sec 

-2ir  rad/sec 

With  these  dynamic  specifications,      the  dynamic  curvature  transform  of 
MECSYN  can  be  used  to  find  the  constraint  links  EF  and  GH  to  constrain 
link  2  relative  to  link  5  in  Fig.  (4-6b).   The  numerical  coordinates  for 
the  chosen  constraints  are 

E  (0.940,  2.313),     F  (0.899,  1.724) 

G  (2.716,  2.475),     H  (1.441,  1.293)  . 
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Step  4.   The  inversion  to  the  Stephenson  2  by  fixing  link  4  from  the 
Stephenson  3  has  all  the  desired  dynamic  specifications  for 
link  2  with  point  D  the  path  generator  and  9  (t)  the  desired 
input  specs.   Figure  (7)  illustrates  the  resulting  Stephenson  2 
mechanism.   The  dashed  path  is  that  for  point  C  of  the  corre- 
sponding inversion  to  the  Stephenson  3  mechanism.   Although  the 
resulting  mechanism  was  satisfying  for  many  reasons,  the  input 
angle  8,  could  not  be  increased  monotonically  even  though  it 
could  be  rotated  completely  in  the  same  closure  for  the  mechan- 
ism.  Therefore,  the  Stephenson  3  mechanism  of  Section  1  synthe- 
sized for  the  same  specifications  is  a  better  mechanism  for  appli- 
cation to  the  stated  transfer  problem.   Note  that  varying  any 
of  the  motion  specifications  a  small  amount  can  dramatically 
affect  the  suitability  of  the  resulting  mechanisms.   Also,  the 
coordinates  U,,V  can  be  changed  arbitrarily  without  changing 
the  basic  transfer  problem.   The  coordinates  of  point  A  were 
varied  many  times  resulting  in  numerous  distinct  Stephenson  2 
mechanisms.   The  final  result  satisfies  more  of  the  necessary 
physical  design  criteria  than  the  others.   The  first  3  steps  of 
this  section  are  illustrated  in  Table  4-1  consisting  of  output 
from  MECSYN. 

As  indicated  in  Chapter  II,  a  Stephenson  1  mechanism  can  also 
be  synthesized  for  the  same  specifications,  being  more  difficult  to 
design  than  a  Stephenson  3,  yet  less  difficult  to  design  than  a 
Stephenson  2.   The  author  concludes  that  for  demanding  specifications 
as  in  these  examples,  a  Stephenson  3  generally  provides  better  prac- 
tical results  for  the  least  effort. 
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Figure  4-7  Transfer  Mechanism  -  Stephenson  2 
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Table  4-1     Stephenson  2   Synthesis 


SVNTH'SIS    OF     MECHANISMS 
FIVfc     MULTIPLE     SfcPAHATEO     POSITIONS 
IN   COPt-AMU    MOTION    -    CASE       3 


0  .  J 
1.0 

SPtCIF 1C*I IUN1 

AI.PHAI0.01         -  2.O00J 

BETAI0.01        '  0.9000 

CAMMAIOtO)        ■  *■«  .0000 

Al.PHA10.il         *  iitSJW 

BETA10.1I         -  .'..-I.:.* 

CA«H«I 0(11  ■  3.1*10 

m.ph*io.2>      -  a.sesB 

BdTA10.il        •  -8.aoS* 

ALPHA!  I. 01         ■  1.9000 

BETA(L.Q)         ■  -1.0000 

6AM«A(I.O)         "  0.0 

Ai.pnAfi.li      ■  uoooo 

BE.TAU.ll       ■  0.0 

GAMMA(I.I)         -  0-0 


1.90O7? 

O.JoIlT 


CUOBOINATSS    UF  THE    CONSTRAINT    LIN*S 
POINTS  ' 

T    »  ItWW  *   »  1.13255 

y    •  2.25209  K    -  0.B7290 


J.BW 


-fc«fc 


-■;0 


SeCONO    CQ6NA 


D3ja.E-t  ocuea 


POS.    DEA. 


2 
3 

57 

MODIFIED    L.L-B 

"iQ          2.59 
2. 55 

L1N«    UN.FH    SAtlO         2. S3 

1MPUT    PABAMETEP.    CONVCBSION 
SIX-L.IN*    AND     ',--:«'.)    <- 
(BASED    ON    STEPHeNSJN-2 

V£-L!N*    MCCHaMSM 

SIA-LINK     -.   ■    .'» HIS*  1 

CASE    NU«BE(*        3 

1 

POINT    A 
2500,           I.2S00) 

POI 
1           0.3730 

T    ,1                                                         POINT    C 

2.07931                 I           0.3*12.          2.29211 

LINK      3 

LENGTH             AN 
0.5*02       161 

3*  JO 

LIN*    • 

l"-i,OI                  ANGLE 

0.9110       U*.*««7 

HE 

A       ? 

STEPHENS  ON -2 
ThETA    5 

THETA     J 

STEPHtNSON-3 
PS  J     3                      PSI    * 

«5 

9 

0 

0000 
1  »16 

0 

a 

• 

0.0 

A. 2832 

0.0 

6. 2* J 2 

1  -1   .3*91 

ft.  21*1 

20.6503 

770.1239 

1.8512 

l  ■■  .  -  1*91              11».»*«T 
I. 1129                -6.2B32 
iJ.S103                  0.0 

0.0                   -155.SS33 
-*. AJ20                -A.2k32 

1 3.2211 

a. 9609 

-2.3609 


-43.9957 

2*6163 


ATCS    OF     THE    CCN 


I".'-'-,     ifl 

1 .19770 


2.4  79  I  9 
-11.992E J 


1.721*0 

0.BA77T 
t.  292-99 


secono  cognate 

-0-914.  1,0731 

1.4*1   ■  1.2931 


DCUBLt-ROCXF-fl 

LIN*     LENC1H     BAI10 
NOG  IF  I  i  n    L.L.B. 


CBANK-eaCKEf* 
LI  Mt  LfNCTM  BATIO 
MOOIFI'O    I   .L.J. 
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c-   Geometric  Specs  Synthesis  of  a  Stephenson  3 

Mechanism  for  Cutting  Large,  Low  Pitch  Gear  Teeth 

This  design  illustrates  the  special  case  of  geometric  synthesis 

for  a  Stephenson  3  six-link  mechanism  as  defined  in  Chapter  II.   The 

specifications  for  the  Stephenson  3  for  dynamic  synthesis  are 

/~00        00        (k)         (k)   ^ 
^Xj(t),      y..(t),      4,  (t),       lh(t)L 

It  was  shown  in  Chapter  II  that  the  dynamic  specs  formulation  contains 
as  a  subset  the  MSP  geometric  synthesis  for  a  link  moving  with  general 
coplanar  motion  relative  to  input  displacement  of  a  crank.   The  subset 
is  obtained  by  replacing  t  with  the  input  *  giving  the  specification  for 
geometric  synthesis  as 


c 


oo         oo         (k>  (kn 

xjW).         y.  (*),         *.(*),  iji    J 


where 

♦  ■  1,      if  =  *'  =  TP".   =  0   . 

The  goal  here  is  to  cut  by  mechanically  controlled  motion  the 
teeth  of  a  very  large,  low  pitch  gear.   The  gear  is  an  8-tooth,  3/8  pitch, 
21.333  inch  pitch  diameter,  6  inch  face  width  involute  spur  gear.   The 
following  quote  comes  from  Ref.  [28]: 

It  is  still  very  difficult  to  cut  a  gear  tooth  coarser  than 
3/4  pitch.   Lufkin  Foundry  and  Machine  Company  have  recently  made 
some  very  well-formed  gear  teeth  twice  this  big  at  3/8  pitch. 

The  gear  teeth  can  be  cut  either  by  oxygen  or  bandsaw  cutting.   With 
oxygen  cutting,  it  may  be  more  efficient  to  use  a  shaped  flame  and 
to  maintain  a  given  flame  orientation  to  the  tooth  surface.   For  band- 
saw  cutting,  it  is  essential  that  the  saw  blade  be  tangent  to  the  tooth 
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surface.   In  either  case,  it  is  desirable  to  have  constant  cutting 
speed  while  the  driving  shaft  angular  velocity  is  constant  to  produce 
a  uniform  cutting  operation.   Also,  a  high  degree  of  accuracy  in  the 
cutting  path,  tangent  angle  and  constant  speed  condition  is  consid- 
ered important. 

A  Stephenson  3  mechanism  will  be  synthesized  to  satisfy  the 
posed  problem.   Higher  order  motion  properties  will  be  derived  from 
the  definition  of  the  input  and  the  involute  curve  and  corresponding 
involute  angle.   The  MSP  specifications  will  be  obtained  for  case  PPP-PP 
to  achieve  close  matching  of  the  involute  geometry.   The  cutting  point 
will  be  point  C  in  the  synthesis  formulation  while  the  angle  <|>  of  the 
corresponding  coupler  BC  will  generate  the  involute  angle  to  maintain 
the  tangency  requirement  for  proper  cutting. 

The  driving  crank  angle  ij,  is  defined  by  the  following  conditions: 

1.  Cutting  speed  is  constant  when  input  speed  is  constant 

dt  *•  r  dt 

where 

r  -  pitch  radius 

K  -  input  scale  factor 

R  -  reference  vector  magnitude  for 
point  on  the  involute 

(f  -  involute  angle. 

2.  Input  conditions  are  initially  zero,  i/in  =  0  at  $  =  0. 
From  condition  1  above 


R 

* 

/ 

dR  ■ 

>    K 

'1 

<ty 

Ro 

*0 
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with  R(0)  =  RQ.   This  shows  that 

R  -  E  =  K  r  t  . 

3.   The  involute  curve  is  described  parametrically  in  terms  of  the 
involute  angle  $  by 

x  =  r(cos  <j>  +  if  sin  $) 
y  =  r(sin  <j>  -   $  cos   (J>) 
such  that 


R  =  /x2+y2  =  r  /H-*2,    R(0)  =  RQ  =  r. 
Hence,  the  input  ij;  is  defined  as 

*  =  £  (R"V  '  1   &**1  "  1)  ■ 

The  first  and  second  derivatives   of   the   functions   <|>(if0,    x(i)i) , 
y(i/0   are  required   for   the  geometric  synthesis   of   case  PPP-PP.      These 
data  are  given  by 


<)>(♦)   -  +  /K2iJ)2  +  2Ki/j  =  g(i)0 
x(\ji)   =  r(cos   g(i)j)   +  g(,j,)    sin  g(i|>)) 
y(W  =  r(sin  g(i(i)   -  gft)   cos  g(«) 
ii  =  *i  fj,1    -  K2»+K 


j|  =  x'dJO    =  rK(K*+l)    cos   g(*),  &  =  y'(,j,)    =  rK(KiM-l)    sin  g0|0 


dip 

d*  [g(*)j3 

—J  =  x"G(,)  =   r{K2    cos  g(^)    -    [K2g(,J,)   +-~]    sin  g(,J,)} 
dl(;z  g\V) 

d2y 

-  y"0|0  =   r{K  sin  gOj,)   +   [K2g(^)   +  -£-]    cos   g^)}    . 

di(/z  gW 
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The  description  of  the  coarse  pitch  involute  spur  gear,  20°  stub 
tooth  by  AGMA  Standards  is: 

Pd  =  21333  =  0-375  =  3/8 

Addendum  -  -^  =  2.133" 
d 

Minimum  =  J_  =  2  667„ 

Dedendum   Pa 
d 

Working  Depth  =  ~-  =   4.267" 
d 

Circular  Tooth  = 
Thickness 

rJ   „  "  («L At.   J  cos  20°  =  10.024" 
base     pitch 

<t>   (at  pitch  circle)  =  20.849° 
P 

4>      (at  addendum  circle)  =  45.497° 

To  obtain  the  scale  factor  K,  choose  the  first  finite  position  j  =  0  (with 
the  PPP)  close  to  the  pitch  point,  say  <j>  *J  21°  and  the  second  finite 
position  j  =  1  (with  the  PP)  close  to  the  addendum  circle,  say  $  sa  40° 
making  A<fQ_>1  =  19°.   Also,  choose  the  crank  rotation  angle  for  these 
positions  to  be  Ai))..  =  75°.   Then  it  follows  that 

*Q  =~   (3.7272°),  f1=j   (12.5813°) 

for  $0   =  21°  and  $  =  40°.   Noting  that  At|f-  .,  ■  75°  the  value  of  K 
becomes  K  =  0.1181.   The  crank  positions  from  the  above  expressions 
using  this  value  of  K  are 

i|>0  =  31.572°,     ,),  =  106.572°  . 
The  synthesis  of  the  mechanism  follows  the  3  steps  outlined  in 
Chapter  II.   These  are: 
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Step  1.   For  the  plane  BC  which  controls  the  cutting  device  the  motion 
specifications  are 


I 

i 

k 

(k) 

X..W 

(k) 
yjW) 

(k) 
♦  .(* 

0 

0 

0 

10.6748 

0.1623 

21.0° 

1 

0 

1 

1.1766 

0.4517 

0.3431 

2 

0 

2 

-0.0245 

0.4513 

-0.2831 

3 

1 

0 

12.1771 

1.0825 

40.0° 

4 

1 

1 

1.1056 

0.9277 

0.2062 

Using  the  dynamic  curvature  transform  in  MECSYN  with  the  replace- 
ment (t  ""   i|/),  synthesize  one  constraint  link  BO   to  control  the 

motion  of  coupler  BC  to  form  the  dyad  CBO.  .   The  numerical  coor- 

b 

dinates  for  the  points  in  this  dyad  are 

C  (10.675,  0.162),   B  (-16.237,  -5.254),   0b  (2.444,  1.204) 

Step  2.   This  step  is  the  solution  of  an  ordinary  geometric  synthesis 

(k)     (k) 
problem.   Specify  the  same  specs  (  x.(ij>)>   y .  (tp) )  as  before 

J       J  (k) 

for  point  C.   The  corresponding  values  for  the  input  specs  <p . 

are  31.572°,  1,  0,  106.572°,  1  for  sequential  positions  I. 

In  the  dynamic  curvature  transform  use  the  replacement  (t  ^  i|») 

to  synthesize  by  MECSYN  two  constraints  on  the  coupler  link  of 

a  four-bar  linkage.   The  cognate  DO  EO  with  coupler  point  C 

00       (k?  e 

has  the  geometric  specs  (  x.(i|i),   y  (\J<) )  at  point  C  coordi- 

(k)   J       J 
nated  with  the  input  i|i.  .   The  numerical  coordinates  of  the 

significant  points  of  this  cognate  are: 
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D  (7.977,  1.731),   0d  (9.208,  1.929) 
E  (10.144,  1.017),  0   (9.552,  4.168) 
where  DO  is  the  Input  crank. 

Step_3.   The  dyad  CBO  of  Step  1  and  the  four-bar  linkage  DO  EO  have 
°  a     e 

the  state  at  point  C  in  common.   Joining  these  systems  at  C 
forms  the  desired  Stephenson  3  gear  cutter  mechanism. 
Figure  (4-8)  illustrates  the  gear  cutter  mechanism  and  the  gear 
geometry.   The  curves  in  Fig.  (4-9)  illustrate  the  quality  of  the  result- 
ing approximation  for  the  x,y  motion  of  the  cutter  along  the  involute 
tooth  and  the  corresponding  error  curves.   The  crank  rotation  angle  iji 
is  used  as  the  abscissa  with  initial  i|i  corresponding  to  $  =  0°  at  the 
base  circle  and  final  *  corresponding  to  <j>  =  46°  on  the  addendum  circle. 
These  curves  show  that  the  involute  is  generated  with  very  high  accur- 
acy except  at  the  ends  of  the  motion.   However,  the  gear  tooth  will  be 
undercut  at  both  ends  consistent  with  the  direction  of  the  error.   If 
the  desired  undercut  is  exceeded,  respacing  the  precision  points  may 
produce  a  better  solution. 

Figure  (4-10)  shows  the  error  curve  for  the  coupler  angle  com- 
pared with  the  true  involute  angle  which  is  surprisingly  small.   It  could 
be  further  improved  by  moving  the  j  =0  location  to  the  left.   The  other 
curve  in  Fig.  (4-10)  shows  the  variation  in  the  velocity  ratio  As/AiJj  which 
should  be  constant  so  that  the  involute  is  cut  at  a  constant  rate.   The 
curve  shows  a  maximum  error  of  ±  1.5%  which  should  prove  very  satisfac- 
tory for  the  intended  purpose. 
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Figure  4-8  Gear  Cutter  Mechanism 
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Figure  4-9  Path  and  Error  of  Gear  Cutter  Mechanism 
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Figure  4-10  Angular  Error  and  Cutting  Rate 


CHAPTER  V 


MECSYN 


If  the  previous  theory  is  to  be  applied  to  the  practical  design 
of  mechanisms,  an  efficient,  complete  computer  program  must  be  available 
due  to  the  lengthy  calculations  involved  and  the  necessity  of  multiple 
solutions  in  obtaining  a  good  design.   This  chapter  presents  a  complete, 
closed  form,  noniterative  computer  program  along  with  a  description  of 
input  and  examples  of  possible  output,  which  performs  all  the  synthesis 
and  analysis  functions  previously  mentioned  and  several  others. 

The  program  is  written  in  Fortran  IV,  avoiding  features  not 
supported  by  ANS  Fortran  so  that  it  is  system  independent  to  a  large 
extent.   Input  and  output  unit  numbers  for  Fortran  are  the  variables 
NI  and  NO  which  are  initialized  in  the  main  program  to  5  and  6. 
The  two  statements  initializing  these  unit  numbers  may  be  changed  to 
reflect  different  input/output  unit  numbers  at  other  installations. 

The  size  of  the  program  when  compiled  on  an  IBM  370/165  with 
Fortran  H  was  about  72K  bytes,  while  the  size  of  the  load  module  created 
was  about  102K  bytes.   Because  of  the  closed  form  nature  of  the  program, 
the  execution  costs  are  nominal.   Execution  of  the  load  module  for  the 
representative  examples  following  the  program  listing  was  about  85  cents 
on  the  University  of  Florida's  IBM  370. 
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As  much  as  possible,  the  output  has  been  left  justified  to 
reduce  output  time  due  to  spacing  on  a  time-sharing  terminal.   Also, 
the  author  has  tried  to  make  the  input  as  simple  as  possible  and  uni- 
form for  all  options. 

In  the  following,  the  terms  "four-position"  and  "five-position" 
refer  to  all  possible  combinations  of  finitely  separated  positions  and 
time  derivatives  up  to  fourth  order  for  those  problems. 

A.   Options  Record 

The  input  for  all  problems  is  preceded  by  an  options  record. 
This  record  indicates  which  options  of  MECSYN  will  be  applied  to  the 
problem  data  in  addition  to  providing  a  title  input  for  each  problem. 
The  format  of  this  record  follows: 

Columns  Description 

0-30      Up  to  six  four-letter  option  words  separated  by 

any  character.   Blanks  and  unrecognized  words  are 
ignored.   The  option  words  may  be  in  any  order. 

31-63  Any  32  character  problem  title. 

66-70  START 

71-75  STOP 

76-80  STEP 

The  last  three  items  are  optional  floating  point  variables  which 
prescribe  the  sweep  through  the  plane  for  solutions  in  a  four-position 
problem.   They  represent  the  initial  angle,  final  angle  and  increment  in 
degrees  (polar  coordinate  angle)  of  the  solution  sweep.   If  they  are 
blank  the  default  values  are  1."  to  171.°  in  steps  of  10.°  which  is 
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equivalent  to  0°-170°  in  the  sweep  (0°il80°).   Since  the  solutions  are 
found  for  each  ray  through  the  origin,  0°-180°  is  sufficient  to  cover 
the  whole  plane.   The  default  values  give  a  maximum  of  36  real  solu- 
tions for  constraint  links. 

There  are  nine  option  words  which  may  he  used  in  various  combi- 
nations.  No  more  than  six  need  to  be  specified  for  any  given  problem. 
If  none  is  specified  the  default  option  is  a  basic  five-position  solu- 
tion.  For  all  other  options  a  five-position  problem  is  solved  unless 
FOUR,  FRFV  or  COGO  is  specified.   The  options  are: 


Option  Words 


COGS 


ANAL 


TRCE 


Option  Description 

With  five-position  solutions,  the  path  cognates, 
Grashof  criteria,  link  length  ratio  are  calcu- 
lated for  all  possible  four-link  mechanisms. 
The  average  length  of  the  two  rays  from  the  pin 
joints  to  the  coupler  point  is  included  In  the 
modified  link  length  ratio  to  give  a  measure  of 
the  distance  from  the  four- link  chain  to  the 
coupler  point. 

When  specified  with  COGS,  this  option  gives  a 
complete  analysis  of  the  mechanisms  created  with 
COGS  for  all  the  specified  precision  conditions. 
Positions  and  derivatives  are  calculated  automat- 
ically for  the  values  of  j  (position)  and  k  (deriva- 
tive) specified  as  input.   This  option  assumes 
the  path  specifications  are  applied  at  the  common 
coupler  point  of  each  ten-link  path  cognate 
mechanism  and  solves  for  the  angular  properties 
of  all  moving  links  in  both  closures  as  described 
in  Chapter  III.   Sometimes  this  information  is 
necessary  to  go  from  one  modular  synthesis  step 
to  another,  in  addition  to  verifying  the  synthesis 
results. 

This  option  traces  the  five-position  solution  and 
writes  intermediate  equations.   Included  are  the 
values  of  the  D   coefficients  used  in  the  problem 

and  the  quartic  polynomial  giving  the  x  coordinates 
of  the  circle  points  for  five  positions. 
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Option  Words 


STLN 


FGEN 


Option  Description 

If  the  path  specifications  satisfy  y  =  ax  +  c  and 
the  path  point  derivatives  satisfy, 

00    (k) 
y At)/   x.(t)  =  a, 

for  five  positions,  a  special  algorithm  must  be 
carried  out  to  obtain  solutions.  This  option 
rearranges  matrices  and  coefficients  to  obtain 
solutions  (using  the  same  basic  equations)  for 
straight  line  motion  specifications.  An  error 
message  notifies  the  user  when  this  option  is 
necessary. 

This  option  converts  four-bar  function  generation 
input  specifications  to  equivalent  specifications 
for  constraint  link  synthesis  and  calls  the  five- 
position  or  four-  and  five-position  constraint 
link  synthesis  subprograms  depending  on  which 
option  is  specified.   This  type  of  function  gener- 
ation does  not  require  inversion  but  angular  or 
path  cognates  as  described  in  Chapter  III. 
If  COGS  and  ANAL  are  chosen  as  options  also,  the 
function  generation  specifications  can  be  observed 
in  the  analysis  of  one  path  cognate  mechanism. 
The  function  generator  can  be  completed  by  angular 
cognates  by  using  the  cognates  only  option  (C0G0) . 
Function  generator  four-bars  are  useful  for  syn- 
thesis problems  for  Stephenson  1  six-bars 
(Refs.  [14,15]).    If  specified  along  with  SIXL, 
this  option  overrides  SIXL. 


SIXL         This  option  performs  input  specification  conver- 
sion similar  to  FGEN  along  with  inversion  to  give 
equivalent  specifications  for  path  related  syn- 
thesis problems  of  geared  five-bars,  Stephenson  2 
six-bars  and  function  generation  for  geared  five- 
bars,  Stephenson  1  and  Stephenson  2  six-bars. 

Function  generation  for  six-link  mechanisms  refers 
to  angular  coordination  of  three  links  for  all 
combinations  of  positions  and  time  derivatives. 
For  path  related  synthesis  problems  a  previous 
constraint  link  synthesis  step  must  be  run  as 
described  in  Chapters  III,  IV  and  Refs.  [14,15 
19].   The  five-position  constraint  link  synthesis 
routine  (or  four-  and  five-position  synthesis 
if  FRFV  is  chosen)  is  called  after  the  SIXL  option 
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Option  Words  Option  Description 

is  completed.   This  option  may  also  be  used  for 
synthesis  of  multi-link  mechanisms  with  more  than 
six  links. 


FOUR 


FRFV 


If  only  a  four-position  problem  is  desired,  this 
option  is  chosen.   The  result  is,  theoretically, 
an  infinity  of  solutions  for  constraint  links 
governed  by  the  sweep  variables  START,  STOP  and 
STEP.   FOUR  overrides  all  other  option  words 
except  COGO.   If  it  is  desired  to  synthesize 
a  six-link  mechanism  for  only  four  positions, 
FRFV  must  be  specified.   That  is,  a  four-position 
and  a  five-position  problem.   Then  options  such 
as  FGEN  and  SIXL  may  be  used  and  solutions  are 
given  for  both  four  and  five  positions. 

This  option  may  be  chosen  along  with  all  other 
options  excluding  FOUR  and  COGO.   The  four- 
position  solution  is  carried  out  first  (after 
any  options  such  as  FGEN  or  SIXL)  using  the  first 
four  sets  of  specifications.   The  five-position 
solution  follows  using  all  five  sets  of  speci- 
fications.  This  implies  one  extra  set  of  fin- 
itely separated  position  specs  or  one  extra  set 
of  velocity  specs  in  addition  to  the  four-position 
specs.   If  FOUR  is  specified  also,  the  next  problem 
will  have  one  extra  input  data  record. 

If  only  path  or  angular  cognates,  Grashof 
criteria  or  link  length  ratios  are  desired,  this 
option  is  specified.   COGO  overrides  all  other 
options.   The  input  consists  of  the  coordinates 
of  five  points  for  any  four-link  mechanism. 
If  angular  cognates  of  a  four-link  mechanism  are 
desired,  the  coupler  point  is  specified  as  one 
of  the  pin  joints.   The  output  will  indicate 
that  cognates,  Grashof  criteria  and  link  length 
ratios  are  desired  but  the  number  of  positions 
specified  will  be  given  as  0.   This  option  allows 
path  cognates  to  be  calculated  separately  for 
selected  four-position  solutions. 

All  the  routines  in  MECSYN  give  solutions  for  any  possible  case 
of  multiply  separated  positions  for  time  dependent  specs.   The  case  num- 
bers are  not  specified  but  are  determined  in  MECSYN  from  the  input  data. 
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There  are  seven  possible  cases  for  five-position  solutions  and  five 
possible  cases  for  four-position  solutions.   These  cases  are  described 
on  comment  cards  in  Subroutine  INPUT.   However,  it  is  not  necessary  to 
be  familiar  with  them  since  they  are  used  only  for  program  control. 
The  parameters  j  (position)  and  k  (derivative)  are  specified  for  each 
set  of  input  data. 

As  many  problems  may  be  run  together  as  desired,  each  preceded 
by  an  options  record  as  long  as  the  correct  number  of  data  cards 
follow  (3,  4,  5  or  6) .   For  each  problem,  a  heading  will  be  printed 
with  "yes"  or  "no"  indicating  which  options  were  selected.   If  options 
are  overriden,  they  will  still  be  indicated  with  a  "yes"  but  the  computa- 
tions will  not  be  carried  out.   The  number  of  specified  positions  will 
be  indicated  by  0,  4,  5  or  45  for  both.   Examples  of  input  and  corre- 
sponding output  follows  the  MECSYN  listing.   Description  of  the  data 
for  each  option  will  be  given  in  the  next  section. 

B.   Input  Data 

All  input  data  following  the  options  record  are  accepted  in  ten- 
column  fields  so  the  same  tab  stops  on  a  time-sharing  terminal  or 
program  on  a  keypunch  may  be  used  for  any  problem.   There  are  three  basic 
types  of  input  depending  on  which  option  is  chosen.   Values  of  the  coun- 
ters j  for  finitely  separated  positions  and  k  for  time  derivatives  begin 
at  zero  and  are  accepted  as  integers  (fixed  point).   All  other  data  are 
accepted  in  real  (floating  point)  format. 

Specifications  are  given  in  a  chosen  fixed  reference  system  and 
all  answers  are  given  in  the  fixed  reference  system  for  the  zero  posi- 
tion of  the  resulting  mechanism.   Derivatives  higher  than  two  must  be 
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specified  for  the  zero  position  (without  loss  of  generality)  but  the 
input  records  containing  the  counters  j  and  k  may  be  placed  in  any  order. 
If  option  FRFV  is  chosen,  however,  the  first  four  j,k  records  are  taken 
for  the  four-position  problem. 

The  description  of  input  records  will  be  given  in  the  order  of 
the  number  required  (3-6)  following  the  options  record. 

If  the  option  COGO  (cognates,  etc.,  only)  is  chosen,  three 
records  follow  the  options  record  giving  the  necessary  information 
about  the  four-link  mechanism. 

Record    Columns  Description 

1        1-10      x  coordinate" 


)Fixed  pivot  of  first  constraint  link 
11-20      y  coordinate, 

21-30      x  coordinate^. 

>  Moving  pivot  of  first  constraint  link 
31-40      y  coordinatej 


1-10      x  coordinate"^ 

Fixed  pivot  of  second  constraint  link 


11-20      y  coordinatej 


21-30      x  coordinate"^ 

>Moving  pivot  of  second  constraint  link 
31-40      y  coordinatej 


1-10      x  coordinate^ 

> Coupler  point  for  path  cognates 
11-20      y  coordinatej 


If  the  option  FOUR  is  chosen  for  a  four-position  problem  only, 
four  records  follow  the  options  record. 
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Record    Columns  Description 

*■"*         1  Counter,  j,  of  finitely  separated  position 

2  Counter,  k,  of  order  of  time  derivative 

3-10  Blank  (to  allow  uniform  tab  stops) 

11-20      x  component^!  Posltion  or  time  derivative  specifi- 
21-30      y  component  >  cations  for  the  coupler  point 
j   (k)     (k) 
(  Xj(t)'   yj(t)) 

31-40      Position  or  time  derivative  specifications  for  an 

(k) 

input,  output  or  coupler  angle  (  9.(t)) 

For  all  other  options  excluding  FGEN  and  SIXL,  five  records  in 
the  same  format  as  those  above  for  option  FOUR,  follow  the  options  record. 

Description 

Finitely  separated  position 

k      k 
Time  derivative,  d  (  )/dt  . 


Record 

Columns 

1-5 

1 

j 

2 

k 

3-10 

Blank 

11-20 

(k) 
x  At) 

21-30 

(k) 
yj(t) 

31-40 

(k) 

e.(t) 

If  FGEN  or  SIXL  is  chosen  the  input  consists  of  the  options 
record,  one  record  giving  the  coordinates  of  three  points,  either  chosen 
arbitrarily  or  solutions  from  a  previous  synthesis  step  and  five  j,k 
records  giving  the  problem  specifications. 
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When  FGEN  is  chosen,  the  next  record  gives  the  coordinates  of 
a  fixed  pivot,  a  moving  pivot  and  a  coupler  point.   These  are  chosen 
arbitrarily  if  only  a  four-bar  problem  is  being  solved.   The  remaining 
records  give  specifications  for  input  and  output  angular  properties. 

Record    Columns  Description 

1        1-10      x  coordinate^ 
11-20      y  coordinate  J 


Fixed  pivot  of  a  constraint  link 


Moving  pivot  of  a  constraint  link 


21-30      x  coordinate 
31-40      y  coordinate 

41-50      x  coordinate  ~\ 

>  Coupler  point 
51-60      y  coordinate  J 

2-6         1        Counter,  j ,  of  finitely  separated  position 

2        Counter,  k,  of  order  of  time  derivative 

3-10      Blank 

11-20      Input  angle  positions  or  time  derivatives, 
(k) 
(  8j(t))   for  the  constraint  link  specified  above 

00 

21-30      Output  angle  positions  or  time  derivatives,  (  $.(t)) 

When  SIXL  is  chosen,  the  next  record  gives  the  coordinates  of 
two  fixed  pivots  and  a  moving  pivot  as  for  a  Stephenson  2  six-bar  or 
a  geared  five-bar  mechanism.   If  a  path  related  synthesis  problem  is 
being  solved,  fixed  pivot  A  may  be  chosen  while  fixed  pivot,  B,  moving 
pivot,  C,  and  one  set  of  angular  specifications  are  given  by  a  previous 
synthesis  step.   If  a  six-bar  function  generation  problem  is  being 
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solved,  points  A,  B  and  C  are  chosen  freely  as  well  as  all  three  sets 
of  angular  specifications. 

A  geared  five-bar  function  generation  problem  allows  choice  of 
A,  B,  C  and  two  sets  of  angular  properties.   The  other  set  is  pre- 
determined by  the  gearing  relationships.   Refer  to  previous  chapters 
for  explanation  of  the  notation  used  here. 

Description 
x  coordinate  of  point  A 
y  coordinate  of  point  A 
x  coordinate  of  point  B 
y  coordinate  of  point  B 
x  coordinate  of  point  C 
y  coordinate  of  point  C 

Counter,  j,  of  finitely  separated  positions 
Counter,  k,  of  order  of  time  derivative 
Blank 

Position  and  time  derivative  specifications  for 
(k) 

92'  9lor  96   (  62j(t)) 


Record 

Column 

1 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

2-6 

1 

2 

3-10 

11-20 

21-30 


31-40 


Position  and  time  derivative  specifications  for  6 
(k)  3 


(k) 

(  83j(t)) 


Position  and  time  derivative  specifications  for  9 
(k)  5 


(k) 

(  85j(t>) 


All  the  above  problems  may  be  solved  for  the  usual  geometric 
synthesis,  also.   Input  for  all  path  properties  will  represent  geometric 
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k       k   k       k 
derivatives  rather  than  time  derivatives  (e.g.,  d  x(t)/dt  -*■  d  x(6)/d9  ). 

Input  for  all  associated  angular  derivatives  will  be  zero  above  first 
order  and  one  for  first  order  (d  6(t)/dt  -+ d  6/d6  ;  6+1;  6, . .  .  ,¥  -+0) . 
For  all  the  problems  involving  angles  only  (FGEN,  SIXL) ,  the  input 
angle  first  derivative  is  specified  as  unity,  while  all  higher  deriva- 
tives are  specified  as  zero.   The  other  angular  properties  represent 
geometric  derivatives  with  respect  to  input  angle  rather  than  time 
derivatives. 

Note  that  all  the  foregoing  synthesis  problems  for  which  time 
dependent  specifications  are  made  must  have  a  specified  time  dependent 
input  which  can  be  obtained  to  drive  the  mechanism,  as  explained  in 
previous  chapters.  For  this  reason,  multi-link  mechanisms  are  best  suited 
to  this  type  of  synthesis  (more  free  parameters).  However,  if  these  con- 
ditions cannot  be  met,  MECSYN  can  still  be  used  for  the  more  restric- 
tive, classical  geometric  synthesis  problems  by  reducing  the  input  as 
shown  above. 

C.   Program  Description 

A  brief  description  of  the  subroutines  in  MECSYN  may  be  useful 
to  the  user.   In  addition  to  this  description,  there  are  several  comment 
statements  within  the  program.   All  real  (floating  point)  variables 
within  the  program  are  of  length  8  (double  precision)  excluding  control 
parameters.   Past  usage  has  shown  this  to  be  necessary  in  many  cases. 
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MAIN 

The  main  program  performs  calls  to  subroutines  but  no  computa- 
tions other  than  rearranging  matrices  for  the  straight  line  problem, 
coordinate  transformations  on  the  constraint  link  solutions  and  deter- 
mining the  center  points  for  the  five-position  solution. 

INPUT 


The  input  data  are  read,  interpreted,  control  paramter  values 
set  and  preliminary  output  is  written.   INPUT  also  calls  subroutines 
according  to  which  options  were  chosen. 

EQUATN 

The  coefficients  of  the  constraint  equations,  D  .,  are  computed 
for  four-  or  five-position  synthesis.   The  coefficients  are  computed 
according  to  the  case  number  which  has  been  calculated  from  the  input 
data  (j  ,k) . 

MATRX 

The  determinant  of  a  (4*  4)  matrix  is  solved  for  the  linearized 
constraint  equations. 

ELIM 

A  variable  is  eliminated  between  two  general  conies  to  obtain 
a  quartic  polynomial  using  Bezout's  method  [23].   The  resulting  quartic 
polynomial  gives  the  x  coordinates  of  the  circle  point  solutions  for  the 
five-position  synthesis  in  the  moving  reference  system.   Coefficients  are 
rearranged  for  the  straight  line  problem  and  the  same  equations  give 
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a  cubic  polynomial  for  the  x  coordinate  of  the  center  point  solutions 
in  the  fixed  reference  system. 

QUARTC 

A  quartic  polynomial  is  rooted  in  closed  form  (real  roots  only)  , 

CUBIC 

A  cubic  polynomial  is  rooted  in  closed  form.   The  resolvent 
cubic  is  solved  for  the  Subroutine  QUARTC. 

COG 

Path  cognates  are  computed  for  a  four-link  mechanism  and  the 
solutions  are  written.   Analysis  and  Grashof  criteria  subroutines  are 
called. 

ANALYS 


All  the  dyads  of  a  ten-link  path  cognates  (Robert's  cognates) 
mechanism  are  analyzed  for  the  number  of  positions  and  time  derivatives 
specified  depending  on  the  case  number.   The  specified  path  properties 
for  synthesis  are  applied  at  the  coupler  point  and  the  angular  proper- 
ties are  found.   This  subroutine  provides  analysis  support  required  for 
synthesis  of  multi-link  mechanisms  and  verifies  the  synthesis  results. 
This  routine  is  called  only  by  COG  (COGS  must  be  chosen  as  an  option). 

INV2 

Inversion  of  a  2 x 2  matrix  is  performed  for  Subroutine  ANALYS. 
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GRASHF 

A  four-link  chain  is  evaluated  according  to  Grashof  criteria. 
Link  length  and  modified  link  length  ratios  are  calculated. 

FOUR 

All  cases  of  four-position  synthesis  are  solved  by  utilizing 
Subroutine  EQUATN.   Four-position  solutions  are  written. 

PET 

The  determinant  of  a  (3<  3)  matrix  is  solved  to  service 
Subroutine  FOUR. 

CRVTRN 

This  routine  calculates  the  fixed  pivot  or  center  point  solu- 
tions for  the  four-position  problem  and  transforms  the  circle  points 
to  the  fixed  reference  system. 

INVERT 


Input  specifications  for  five-  and  six-link  synthesis  and  four- 
link  function  generation  are  converted  to  specifications  for  constraint 
link  synthesis  by  inversion  and  analysis.   The  four-  or  five-position 
synthesis  routines  are  called  by  using  these  specifications. 

A  listing  of  MECSYN  along  with  examples  of  input  and  output 
follows . 
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CHAPTER  VI 


RECOMMENDATIONS  AND  CONCLUSIONS 


There  are  several  subjects  which  may  be  logical  extensions  of 
the  preceding  study.   Among  these  are  the  time  dependent  analogy  of 
classical  curvature  theory  and  design  optimization  techniques  which  may 
be  applied  to  the  synthesis  of  multi-link  mechanisms.   Before  concluding, 
preliminary  formulation  of  the  first  topic  will  be  given  along  with 
recommendations  for  further  research  on  both  topics. 

A.   Time  Dependent  Analogy 

of  Classical  Curvature  Theory 

Classical  curvature  theory  is  concerned  with  the  specification 
of  the  instantaneous  motion  of  a  planar  body  for  geometric  properties 
such  as  tangent,  curvature,  rate  of  change  of  curvature  and  the  deter- 
mination from  this  information  several  special  loci  within  the  plane 
that  may  be  used  for  application  of  constraints.   This  type  of  synthesis 
is  known  as  infinitesimally  separated  position  (or  ISP)  synthesis. 
Some  of  these  special  loci  are,  for  example,  the  locus  of  all  points 
which  have  zero  curvature  (3  ISP)  so  that  a  slider  or  prismatic  pair 
constraint  may  be  used.   In  addition,  the  locus  of  all  points  having 
stationary  curvature  (with  respect  to  the  geometric  parameter  9  for 
4  ISP)  can  be  obtained  so  that  a  revolute  pair  or  pin  jointed  constraint 
link  may  be  used. 
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The  most  definitive  works  on  this  subject  have  assumed  8  ~   t   so 
all  time  derivatives  are  eliminated  (6  =  1,  6  =  0,  6  =  G,  "&'=  0)  in  addi- 
tion to  the  assumption  of  a  special  reference  x  =  y  =  9  =  0  and  x  =  y  =  0. 
The  second  condition  means  the  first  position  is  at  the  velocity  pole. 
Also,  in  classical  curvature  theory  all  derivatives  are  geometric  deriv- 
atives  x'=dx/d6,  etc.   Here,  all  derivatives  will  be  with  respect  to 
time  including  the  derivatives  of  8  (x=dx/dt,  9 = d29/dt  ,  etc.). 
To  allow  this  formulation  to  be  simplified  to  a  reasonable  extent,  the 
same  assumption  will  be  made  by  placing  the  first  position  at  the 
origin  (0  =  x  =  y  =  O)  but  the  first  position  will  not  be  coincident  with 
the  velocity  pole  (x  j   0,  y  j   0)  and  (6  4   0,  6  +   0,  9'  4   0,  T 4   0)  will 
be  considered.   The  analogous  loci  to  some  classical  loci  from  curva- 
ture theory  will  be  shown  in  addition  to  loci  not  found  in  this  geometric 
based  theory. 

As  before,  the  specifications  are  given  by 

k    ,k    k   ,k     k    k 

d  x      d  y     n    d  9 

X  =  ~~k  •  y=   k  "   9  =  ~ k  * 
dt*      dtk        dt 

The  fixed  reference  system  is  given  by  (U,V)  while  the  reference  system 
fixed  in  the  moving  body  is  (u,v).  With  the  conditions  x  =  y  =  9  =  0, 
the  coordinate  transformation  and  time  derivatives  are  given  by: 
U  =  u  V  =  v 

ti  =  -9v  +  x  V  =  9u  +  y 

(6-1) 

U  =  -02u  -  6*v  +  x  V  =  9u  -  02v  +  y 

U'  =  -39  9u  +  93v  -  "e  v  +  x  V  =  -03u  +  *9*  u  -  39  9v  +  y 
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k  k 
The  derivatives  of  the  path  specifications  (x,y)  may  be  expressed 

in  terms  of  the  equivalent  geometric  derivatives  and  time  derivatives 

of  theta  by  the  chain  rule  for  derivatives  to  facilitate  comparison 

with  classical  curvature  theory  loci. 

Since  this  is  a  discussion  of  possible  further  study,  the 

properties  will  be  extended  to  the  third  order.   Fourth  order  loci  may 

be  obtained  in  the  same  fashion.   Also,  for  these  special  properties 

or  loci  to  be  useful  to  the  designer  using  MECSYN,  they  must  be  written 

in  terms  of  the  coefficients  D  .  given  in  Chapter  II  so  that  they  are 

mx. 

completely  general  (x  4   0,  y  4   0,  6  ^  0) . 

Some  preliminary  formulas  will  prove  valuable  for  finding  the 

desired  loci.   After  the  specifications  for  the  origin  of  the  (u,v) 

k  k  k 
system  are  given,  (x,y,6),  the  path  of  any  point  (u,v)  at  the  point 


(U,V)  has  the  curvature 

_  [by  -  iiv] 


(6-2) 


[i2  ♦  V2]3/2    " 

The  acceleration  of  a  point  (u,v)  at  the  point  (U,V)  moving  along 
a  path  S  is 

£(t)  -■£&!+   ff]2N.  (6-3) 

dt2  "    \dtl     ~ 

where  T_  is  the  unit  vector  tangent  to  the  path  and  N  is  the  unit 
vector  normal  to  the  path  at  (U,V)  where 

(df f  =  U2  +  V2  •  (6-4) 
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The  jerk  of  a  point  (u,v)  at  (U,V)  is  given  by 


i(t)  -£«<t> 


^ 


T  + 


Ids \.  dK  dS   d^S 

UtJ    dt       JK  dt       2 


(6-5) 


The  fourth  derivative  (jounce)  may  be  used  to  find  the  tangen- 
tial and  normal  loci  in  the  same  manner  as  for  the  second  and  third 
order  properties.   The  k   order  motion  poles  may  be  of  interest  to  the 
designer.   They  are  defined  as  the  points  in  the  plane,  (u,  ,v,),  which 
have  no  k   order  motion  vector  (e.g.,  zero  velocity,  zero  acceler- 
ation, etc.) 


I  =  /  U  + 


■*     U  =  0,   V  =  0 


(6-6) 


The  two  poles  P  and  P  will  be  shown  to  illustrate  the  dependence  of 

k 
the  higher  poles  on  6 . 

Using  the  chain  rule  for  derivatives  where  g  5  9/§2, 
P.      U.  =  -y/6  =  -y' 


vx  =  x/e 


(6-7) 


P2      "2  =  (-ye+xS^/d^+B2)  =  (-y"g    -y'g2+x"+x'g)/(l+g2) 

v2  =  (e2y+exy(e1,+e2)   =  (y"+y'g+x"g  +x-g2)/(1+g2)  . 

Note  that  P  is  dependent  only  on  the  geometric  derivatives,  while 

P-  is  dependent  also  on  9, 

curvature  theory  set  6  =  1,  9  =  0  which  gives 


For  the  analogy  of  P  in  classical 
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Cubic  of  Stationary  Curvature 

A  locus  of  considerable  value  in  geometric  ISP  synthesis  is  the 
locus  of  stationary  curvature  (4  ISP).   Although  MECSYN  contains  rou- 
tines for  solutions  for  all  combinations  of  4  FSP  and  4  ISP  (4  MSP)  for 
both  time  dependent  and  geometric  specifications,  this  locus  will  be 
given  here  to  show  its  relation  to  the  other  special  loci  under  con- 
sideration.  The  problem  can  be  stated  as: 

Given  a  set  of  specifications  for  the  third  order  motion  of 

kkk 
a  body  (x,y,6),  find  the  locus  of  all  points  (u,v)  which  have  at  the 

points  (U,V),  stationary  curvature  (with  respect  to  time). 

Since  the  condition  of  stationary  curvature  means  the  points 
(u,v)  have  third  order  contact  with  a  circle,  these  points  may  be  used 
as  pin  joints  of  circular  constraint  links  for  the  motion  as  specified. 
In  classical  curvature  theory,  this  locus  is  expressed  in  terms  of 
geometric  derivatives  only.   It  will  be  shown  (as  in  Chapter  II)  that 
this  locus  is  not  sufficient  to  synthesize  a  mechanism  for  third  order 
time  derivatives.   The  additional  condition  will  be  shown  to  be  related 
to  the  tangential  jerk  component.   The  constraint  equation   formulation 
in  Chapter  II  shows  the  necessity  of  this  extra  condition  for  all  cases 
of  time  dependent  MSP. 

The  locus  of  stationary  curvature  is  obtained  by  setting  the 
time  derivative  of  Eq.  (6-2)  equal  to  zero. 

d<  _  [u2+v2]  Ei'v-u  id  -  3  [fiv-uv]  [uu+vv] 

It "   °  (6-8) 

[A*  +  V2]5/2 
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or 

[u2+v2]  [u  v  -  'if  v]  -  3[uv  -  ttvj  [uti  +  irf\  -  a  (6-9) 

if  U2  +  V2  ?«  O  .  (6-10) 

Equation  (6-9)  appears  to  be  fourth  order  in  (u,v)  but  the  analogy  in 
classical  curvature  theory  has  shown  this  to  be  cubic  and  is  called 
the  cubic  of  stationary  curvature.   The  condition,  inequality  (6-10), 
will  be  necessary  for  some  of  the  following  loci.   With  the  usual 
assumption  that  the  origin  of  (u,v)  is  placed  at  P.,  this  condition 
reduces  to  u  ^  0,  v  4   0.   However,  it  represents,  in  general,  the 
exclusion  of  all  points  (u,v)  which  lie  at  the  intersection  of  the 
locus  in  question  and  a  circle  represented  by  inequality  (6-10); 

(u2+v2)  +  2(y/9)u  -  2(x/6)v  +  (x/6)2  +  (y/6)2  t   0 
or 

(u2+v2)  +  2y'u  -  2x'v  +  x'2  +  y'2  ^  0  .  (6-11) 

For  all  the  loci  to  which  this  condition  pertains,  one  solution  is 
the  pole  P  =  (-y',x'). 

Inflection  Circle 


The  locus  of  all  points  (u,v)  having  a  characteristic 
curvature,  k,  is 

[uV  -  UV]  -  <[u2  +  V2]3/2  =  0  .  (6-12) 

This  is  sixth  order,  in  general.   For  k  =   0  (useful  for  3  ISP  sliding 
constraint)  then 

UV  -  UV  =  0    if  U2  +  V2  ^  0  .  (6-13) 
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Note  that  this  is  given  by  the  third  term  of  the  cubic  of  stationary 
curvature  and  is  referred  to  as  the  inflection  circle.   Any  point  (u,v) 
on  this  circle  undergoes  an  inflection  for  the  given  motion  specifi- 
cations.  Using  the  chain  rule,  the  inflection  circle  can  be  written  in 
terms  of  geometric  derivatives  only,  as 

u2+v2  +  (y'-x")u  +  (x'-y")v  +  x'y"  -  x'V  =  0  (6-14) 

with  center  at 

[(-y'+x")/2,  (-x'+y")/2] 
and  radius 

[x"y'-x'y"  +  (y'-x")2/4  +  (x'-y")2M]%  . 

Zero  Tangential  Acceleration  Circle 

From  Eq.  (6-3),  the  tangential  acceleration  is 

at  =  d2S/dt2  . 

The  locus  of  all  points  (u,v)  having  a  characteristic  tangential 
acceleration,  a  ,  is 

[fit)  +  VV]  -  a  [U2  +  V2]2  =  0  .  (6-15) 

This  is  fourth  order,  in  general.   For  a  =0,  the  zero  tangential 
acceleration  circle  arises. 

uii  +  W  =  0   if  U2  +  V2  )i  0.  (6-16) 

Note  that  this  is  given  by  the  fourth  term  of  the  cubic  of  stationary 
curvature.   Equation  (6-16)  may  be  written 
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g(u2+v2)  +  (y"+2gy'-x')u  -  (x"+2gx'+y ' )v 

+  x'x"  +  y'y"  +  g(x'2+y'2)  =  0  (6-17) 

where 

g  2  9/82  . 

This  shows  that  Eq.  (6-16)  is  dependent  on  the  time  functions  8  and  6. 
The  analogy  to  Eq.  (6-17)  in  classical  curvature  theory  may  be  found 
by  setting  6  =  1,  8  =  0,  giving 

(y"-x')u  -  (x"+y')v  +  x'x"  +  y'y"  =  0  .  (6-18) 

This  is  the  equation  of  the  normal  to  the  inflection  circle  passing 
through  the  pole  P_  =  (-y',x'). 

Zero  Normal  Acceleration  Circle 

The  normal  acceleration  is  given  by  Eq.  (6-3)  as 
,  2 


an  =  K(f) 


The  locus  of  all  points  (u,v)  having  a  characteristic  normal  acceler- 
ation, a  ,  is 
'   n' 

[UV  -  UV]  -  an[U2  +  V2]%  =  0   .  (6-19) 

For  a  =  0  this  locus  is  identical  to  the  inflection  circle.   However, 
for  a  jl  0  the  two  loci  are  distinct.   Recall  the  characteristic 
curvature  locus  is 

[UV  -  UV]  -  k[U2  +  V2]3/2  =  0 
which  is  sixth  order,  while  Eq.  (6-19)  is  at  most  fourth  order. 
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Zero  Normal  Jerk  Circle 

The  normal  jerk  is  given  by  Eq.  (6-5)  as 

,  2 

d<  ,  ,  dS 

- f  3k  — 

At-  At- 

dt2 

but 


4  =  (M\     iE.  +  •»  dS  d£s 
]n   Idtj   dt      dt   2 

fdS^2  dan   dK  /-dS^j 


2  da     ,  /  \  2  , 

dt  ,4 


dt 
So  j  may  be  written 

■   =  d3n  ,   dS  d2S 
-"n   dt   K  dt  ,2   ' 

From  the  above  it  can  be  shown  that  the  locus  of  all  points  (u,v) 

which  have  a  characteristic  normal  ierk,  i  ,  is 

j    .  Jn> 

[U  V  -  U  V]  -  Jn[U2  +  V2]^  =  0  .  (6-20) 

For  j   =  0,  the  zero  normal  jerk  circle  arises 

U  V  -  U  V  =  0    if  U2  +  V2  jf  0  .  (6-21) 

This  equation  is  represented  by  the  second  term  of  the  cubic  of 
stationary  curvature. 

Zero  Tangential  Jerk  Circle 

The  tangential  jerk  is  given  by  Eq.  (6-5)  as, 
d3s 


.3c     '-^ 


^         dt3 
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The  locus  of  all  points  (u,v)  which  have  a  characteristic  tangential 
jerk,  j  ,  is 

[U  V   +  V  V]  -  Jt[U2  +  V2]^  =  0  .  (6-22) 

When  Jt  "  0.  Eq.  (6-22)  gives  the  zero  tangential  jerk  circle. 

U  U  +  V  V  =  0    if  U2  +  V2  4   0  .  (6-23) 

Note  that  this  equation  is  not  represented  in  the  cubic  of  stationary 
curvature.   For  a  geometric  problem,  it  would  not  be  necessary,  but 
for  a  problem  involving  specifications  which  are  a  function  of  time 
this  equation  represents  the  essential  tangential  properties  of  the 
constraint  circle  point  moving  along  the  path.   Thus  the  second  condi- 
tion for  synthesis  on  page  29  is  the  same  as  the  inclusion  of  this 
condition.   Both  signify  the  need  for  a  synthesized  input  which  effec- 
tively guarantees  the  constraint  point  or  circle  point  will  have  the 
necessary  tangential  motion  properties  needed  to  satisfy  the 
specifications . 

Other  loci  such  as  the  characteristic  normal  and  tangential 
jounce  loci  may  be  found  by  differentiating  Eq.  (6-5)  with  respect  to 
time.   An  important  point  to  the  designer  is  Ball's  point.   In  clas- 
sical curvature  theory,  Ball's  point  is  a  point  where  the  curvature 
of  the  path  has  the  stationary  value  zero  or  third  order  contact  with 
a  straight  line.   This  point  may  be  used  for  a  4  ISP  linear  or  sliding 
constraint  similar  to  the  3  ISP  sliding  constraints  of  the  inflection 
circle.   The  analogy  of  this  point  may  be  found  by  the  intersection  of 
the  zero  normal  acceleration  circle  (inflection  circle)  and  the  zero 
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normal  jerk  circle.   Also  it  may  be  of  interest  to  investigate  the 
meaning  of  the  intersection  of  other  loci  of  characteristic  acceler- 
ations and  jerks. 

The  characteristic  acceleration  and  jerk  loci  appear  to  be 
fourth  order,  while  the  zero  acceleration  and  jerk  loci  are  all  second 
order.   It  may  be  of  interest  to  investigate  the  nature  of  these 
characteristic  loci  and  their  possible  use  to  aid  the  designer.   The 
zero  loci  pertain  to  two  common  types  of  constraints,  i.e.,  circular 
and  linear,  while  the  characteristic  loci  may  be  used  to  describe  any 
type  of  general  constraint. 

B.   Mechanism  Evaluation  for 
Design  Optimization 

The  term  optimization  used  in  the  context  of  mechanism  design 
does  not  refer  to  classical  optimization  methods  such  as  those  based 
on  calculus  of  variations  but  rather  to  the  determination  of  effective, 
easily  computed  criteria  by  which  synthesized  mechanisms  may  be  rejected 
or  accepted.   There  are  exceptions  to  this,  such  as  optimal  precision 
point  spacing  for  function  generation  using  Chebyshev  polynomials, 
a  method  well  known  to  kinematicians.   However,  this  is  an  optimization 
method  that  is  applied  before  the  synthesis  of  the  mechanism. 

Although  the  synthesis  process  gives  the  designer  great  freedom 
in  the  choice  of  motion  specifications,  frequently  the  resulting  mechan- 
ism must  also  satisfy  other  conditions  such  as  force  transmission, 
mobility  and  size  requirements.   Some  of  these  criteria  are  time  consum- 
ing to  calculate,  especially  when  many  mechanism  solutions  can  be  gen- 
erated for  the  same  specifications  by  variation  of  a  parameter  (as  with 
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most  six-link  mechanisms) ,  or  by  relaxing  the  specifications  (four- 
position  problem).   Accordingly,  some  of  these  criteria  have  been 
included  in  MECSYN.   Other  criteria  are  not  well  developed,  e.g., 
general  mobility  criteria  for  six-link  mechanisms  similar  to  those 
developed  by  Grashof  for  a  four-link  chain.   Some  of  these  criteria 
may  be  determined  by  a  complete  motion  analysis  of  the  synthesized 
mechanisms,  a  very  undesirable  procedure  considering  computer  time  and 
necessary  core  storage. 

After  all  the  necessary  criteria  have  been  formulated  and 
included  in  MECSYN,  the  designer  may  wish  to  compare  these  criteria 
personally  or  in  the  event  of  multiple  solutions  a  cost  function  or 
penalty  function  may  be  used  to  automatically  eliminate  mechanisms. 
However,  most  mechanism  applications  require  different  criteria  for 
design.   Thus,  any  reasonable  scheme  for  a  penalty  function  must  allow 
the  designer  to  input  weight  factors  for  each  criterion.   A  short 
description  of  several  design  criteria  of  importance  in  multi-link 
synthesis  follows. 

A  most  obvious  condition  for  design  is  the  relative  size  of 
a  mechanism.   Computation  of  the  link  lengths  and  link  length  ratios 
(ratio  of  longest  to  shortest  link)  tells  the  designer  if  the  mechan- 
ism fits  in  the  allowable  space  and  if  the  mechanism  can  be  effectively 
constructed  (very  small  links  may  create  manufacturing  problems) . 
MECSYN  computes  link  lengths,  link  length  ratios  and  modified  link 
length  ratios  for  all  four-link  modules  and  their  cognates.   This  may 
easily  be  extended  to  compute  link  length  ratios  for  a  complete  multi- 
link  mechanism.   The  modified  link  length  ratio  uses  an  equivalent 
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link  length  for  a  ternary  coupler  link.  The  average  of  the  two  rays 
extending  from  the  pin  joints  to  the  coupler  point  is  used  as  an  addi- 
tional link  in  the  link  length  ratio  calculation. 

A  useful  measure  of  mobility  of  a  four-bar  mechanism  is 
Grashof's  criteria  which  give  several  rules  for  the  determination  of 
how  the  input  and  output  cranks  rotate  from  the  mechanism  proportions. 
This  allows  immediate  evaluation  of  whether  or  not  a  mechanism  may 
have  a  continuously  rotating  input  crank,  for  example.   These  criteria 
are  evaluated  for  all  four-bar  mechanisms  and  their  cognates  in  MECSYN. 
At  present,  this  can  be  accomplished  for  six-link  mechanisms  only  by 
a  complete  motion  analysis.   Criteria  similar  to  Grashof's  for  six- 
link  mechanisms  would  be  the  most  desirable  way  of  providing  this  mobil- 
ity information.   An  additional  problem  arises  in  six-link  mechanisms 
for  which  the  input  is  a  ternary  link.   Frequently,  although  the  input 
may  rotate  full  cycle  in  the  course  of  the  mechanism's  motion,  this 
input  rotation  is  not  monotonic  increasing  but  requires  several  rever- 
sals.   Since  this  would  not  allow  practical  implementation  of  the 
mechanism,  good  mobility  criteria  should  detect  this  condition. 

A  very  important  criterion  for  evaluation  of  mechanisms  is  the 
force  transmission  property  of  the  mechanism;  i.e.,  what  portion  of  the 
force  at  the  input  is  transmitted  to  the  output.   In  a  four-bar  mechan- 
ism, this  property  is  clearly  defined  by  the  transmission  angle  (angle 
between  the  coupler  and  output  link) .   For  some  multi-link  mechanisms 
it  is  more  difficult  to  find,  and,  in  general,  can  be  determined  by 
calculating  the  influence  coefficients  or  velocity  ratios  between  the 
input  and  output  links. 
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Although  the  two  or  three  modules  comprising  a  five-  or  six- 
link  mechanism  are  guaranteed  to  be  compatible  for  the  precision 
points  specified  by  the  synthesis  process,  occasionally  they  are  not 
compatible  for  the  full  motion  of  one  module.   Again,  this  could  be 
determined  by  a  coupler  curve  analysis  of  one  portion  of  the  mechanism 
and  a  point-by-point  check  on  the  remaining  portions  of  the  mechanism 
to  assure  they  will  extend  to  these  limits.   For  many  mechanisms,  as 
may  arise  from  a  four-position  problem,  this  could  be  expensive  in 
terms  of  computer  time.   This  might  be  avoided  by  an  efficient  method 
of  determining  coupler  curve  bounds  without  actually  computing  the 
coupler  curve.   Then  dyads,  for  example,  may  be  checked  against  these 
bounds  to  determine  if  they  are  compatible  for  the  complete  cycle. 

A  similar  problem  arises  when  a  dyad  needs  to  be  in  different 
closures  for  different,  finitely  separated  positions.   The  MECSYN  anal- 
ysis of  dyads  at  precision  conditions  gives  the  closure  a  dyad  must  be 
in  to  satisfy  the  specifications.   Determination  of  the  dyad  closures 
for  a  series  of  points  on  a  coupler  curve  between  the  specified  fin- 
itely separated  positions  will  indicate  whether  or  not  the  dyad  is  in 
the  correct  closures.   If  a  four-bar  module  of  a  multi-link  mechanism 
satisfies  the  specifications  in  two  different  closures,  it  can  be 
iramediatly  discarded.   This  type  of  closure  test  would  be  useful  but 
is  not  provided  by  MECSYN  at  present. 

Another  similar  problem  is  that  of  the  order  in  which  finitely 
separated  position  specifications  are  satisfied.   This  is  not  a  fre- 
quent problem,  if  the  specifications  are  chosen  to  give  a  smooth  pro- 
gression of  finitely  separated  position. 
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The  preceding  two  design  difficulties  are  less  likely  to  arise 
when  higher  order  properties  are  specified.  They  occur  mainly  when 
many  finitely  separated  positions  are  specified. 

When  a  function  or  functions  are  approximated  by  several 
precision  conditions  as  in  the  gear  cutter  example  of  Chapter  IV, 
it  is  important  to  know  the  error,  that  is,  the  difference  between  the 
desired  functions  and  the  generated  functions.   This  requires  an  anal- 
ysis of  the  synthesized  mechanism  for  the  desired  functions. 

By  the  methods  presented  here,  the  complete  motion  of  a  body 
may  be  coordinated  with  an  input.   Thus  an  angular  function  and  a  path 
function  are  both  generated.   The  well-known  method  of  optimal  preci- 
sion point  spacing  by  Chebyshev  polynomials  should  have  an  analogy 
for  optimal  precision  position  respacing  for  synthesis  of  path  and 
angular  function  generators.   A  theory  such  as  this  would  be  very 
useful  in  problems  such  as  the  gear  cutter  mechanism  of  Chapter  IV. 

If  all  the  preceding  suggestions  are  implemented  in  MECSYN, 
the  designer's  task  would  be  gr-jatly  reduced,  leaving  only  the  problem 
of  choosing  the  desired  specifications  and  design  criteria. 

C.   Conclusion 

The  theory  and  design  methods  presented  herein  allow  the 
synthesis  of  a  large  class  of  mechanisms  for  more  realistic  problem 
specifications  than  available  with  geometric  synthesis  techniques. 
The  author's  previous  work  in  five-  and  six-link  mechanism  synthesis 
methods  has  been  extended  to  include  many  previously  unsolved  geometric 
synthesis  problems. 
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All  the  synthesis  problems  and  all  cases  of  mixed  specifica- 
tions presented  can  be  solved  by  the  computer  program  MESCYN  in  addi- 
tion to  other  problems  not  treated  here. 

MECSYN  is  a  completely  closed  form,  noniterative,  inexpensive 
program  that  may  be  used  on  most  computers.   MECSYN  can  be  used  effec- 
tively with  minor  knowledge  of  the  underlying  theories. 

Hopefully,  the  methods  and  computational  tool  presented  here 
will  encourage  the  solution  of  practical  problems  by  mechanism  synthesis. 
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